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ABSTRACT 

The specific frequencies of globular cluster systems, Sn oc A/tot/iv,gai oc Mgcs/Afstars, are discussed in 
On , terms of their connection to the efficiency of globular cluster formation in galaxy halos, which is claimed 

On ■ to reflect a generic aspect of the star formation process as it operates even at the current epoch. It is 

On ' demonstrated that the total masses of GCSs are little affected by the dynamical destruction of low-mass 

clusters at small galactocentric radii. This permits direct, empirical estimates of the cluster formation 
■ efficiency by mass, ed = M^]!^ / , even after 10^° years of CCS evolution. However, the standard 

I practice of using only the stellar luminosities of galaxies as indicators of their initial total gas masses 

(and thus relating Sn to ed in one step) leads to serious conceptual problems, that are reviewed here. 
The first specific frequency problem, which is the well known tendency for many brightest cluster galaxies 
O^l i to have higher than average Sn, is a global one; the second specific frequency problem is a local one, 

in which the more extended spatial distribution of GCSs relative to halo stars in some (not all) bright 
ellipticals leads to Sn values that increase with radius inside the galaxies. Extending similar suggestions 
in the recent literature, it is argued that these trends in Sn do not reflect any such behavior in the 
underlying ed; rather, both of these problems stem from neglecting the hot. X-ray emitting gas in and 
' around many large ellipticals, and both may be alleviated by including this component in estimates of 

' This claim is checked and confirmed in each of M87, M49, and NGC 1399, all of which have been 

On ! thought to suffer from one or the other of these Sn problems. Existing data are combined to construct 

On ■ GCS surface density profiles that extend over nearly the whole extents of these three galaxies, and a non- 

parametric, geometrical deprojection algorithm is developed to afford a direct comparison between the 
volume density profiles of their GCSs, stars, and gas. It is found, in each case, that pd oc (pgas + Pstars) at 
radii beyond roughly a stellar effective radius, inside of which dynamical evolution may have depeleted 
O I the initial GCSs. The constant of proportionality is the same in all three galaxies: ed = 0.0026 ± 0.0005. 

Taken together, these results suggest that GCSs generally should be more spatially extended than stellar 
^ ' halos only in gas-rich galaxies that also have a high global specific frequency. 

The implication that ed might have had a universal value is supported by global GCS data for a sample 
of 97 giant ellipticals, brightest cluster galaxies, and faint dwarfs. The total globular cluster populations 
in all of these early-type systems are in excellent agreement with the predictions of a constant ed at the 
, level observed directly in M87, M49, and NGC 1399; all systematic variations in GCS specific frequency 

' between galaxies are shown to result entirely from different relations, in different magnitude ranges, 

between M™* and the present-day Lv,gai- An identical ed is also calculated for the Pop. II spheroid of 
the Milky Way, and is indicated (although less conclusively) for the ongoing formation of open clusters. 
The inferred universal cluster formation efficiency, of ~0.25% by mass, should serve as a strong constraint 
on general theories of star and cluster formation. The associated inference of a non-universal formation 
efficiency for unclustered stars is considered, particularly in terms of the suggestion that this might 
result, both in dwarf galaxies and at large galactocentric radii in the brightest ellipticals, from feedback 
and galactic winds. Implications for a merger-formation model of early-type GCSs, and for the proposed 
existence of intergalactic globulars in clusters of galaxies, are briefly discussed. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: individual (M87, M49, NGC 1399) 
— galaxies: star clusters — globular clusters: general — stars: formation 

1. INTRODUCTION 

As the study of galaxy formation comes more and more to require a reliable description of star-gas interactions and 
feedback on relatively small spatial scales, it is forced towards questions that overlap with fundamental issues in star 
formation theory. Conversely, as evidence accumulates that most stars today — in environments ranging from the disk of 
the Milky Way to circumnuclear starbursts — are born not in isolation but in groups, the study of star formation in general 
seems increasingly to demand a focus on larger scales than might previously have been considered. Although these two 
disciplines have not yet converged — the smallest scales currently of importance for galaxy formation can still be one or 
two orders of magnitude greater than the "large" scales in star formation — they will necessarily continue to approach one 
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another as further progress is made. It is therefore important to identify and understand those areas of truly common 
interest, on the largest stellar scales and the smallest galactic ones. In many ways, the systems of old globular clusters in 
galaxy halos are positioned just at this interface; they can potentially be exploited to constrain aspects of both local and 
global star-formation processes, in both protogalactic and present-day settings. 

It has, indeed, long been recognized that the global properties of globular cluster systems (GCSs) might be used to good 
effect as probes of galaxy formation and evolution. (This is especially true of bright early-type galaxies, whose GCSs are 
more populous and better studied than those in typical spirals like the Milky Way.) For example, following the suggestion 
that large fractions of the GCSs in ellipticals may have formed in major merger events (e.g., Schweizer 1987; Ashman & 
Zepf 1992; Zepf fc Ashman 1993), and the related discovery of young super star clusters in systems like the Antennae 
galaxies (Whitmore & Schweizer 1995), much recent discussion has centered on the interpretation of the distribution of 
metallicities (or broad-band colors) among E-galaxy globulars as clues to their hosts' dynamical histories (see Geisler, 
Lee, & Kim 1996; Forbes, Brodie & Grillmair 1997; Cote, Marzkc, & West 1998; Kissler-Patig, Forbes, & Minniti 1998). 
A somewhat older approach has used just the total population of GCSs, as a function of parent galaxy luminosity, Hubble 
type, and local environment, to make arguments relating to issues such as biasing in a cold dark matter cosmology (West 
1993); the merger history of ellipticals (Harris 1981; van den Bergh 1984; Ashman & Zepf 1992; Forbes et al. 1997); 
the formation of cD galaxies and their envelopes (McLaughlin, Harris, & Hanes 1993, 1994); star formation in galaxies 
at the centers of cluster cooling flows (Harris, Pritchet, & McClure 1995); the influence of early galactic winds on the 
evolution of giant ellipticals (Harris, Harris, & McLaughlin 1998); and even larger-scale aspects of galaxy clusters such 
as populations of intergalactic stars (e.g.. West et al. 1995; Harris et al. 1998). Yet a third characteristic of GCSs — their 
radial distributions (the projected number density of globulars around a galaxy as a function of galacto centric radius) has 
also served as the basis for discussions of elliptical galaxy formation (e.g., Harris 1986; Ashman & Zepf 1992; Kissler-Patig 
1997; van den Bergh 1998). 

Rather less developed, however, is a clear sense of the relevance of GCSs to local star formation on smaller scales: What 
relation, if any, do the fundamental properties of globular cluster systems bear to the way stars form now, in the dense 
clumps within Galactic giant molecular clouds? Can the study of one of these phenomena possibly shed light on the other? 
Larson (1988, 1993, 1996), has emphasized the possible parallels between the basic pattern of local star formation and the 
birth of globular clusters within large protogalactic clouds in a scenario for the collapse of galaxies from clumpy initial 
conditions (as in, e.g., the classic picture of Searle & Zinn 1978, and current models of hierarchical galaxy formation). 
A more quantitative discussion along these lines, including an exposition of links to several of the gross characteristics 
of GCSs, is given by Harris & Pudritz (1994). Most recently, McLaughlin & Pudritz (1996) and Elmegreen & Efremov 
(1997) have focussed in detail on the globular cluster luminosity function (the number of clusters per unit magnitude, 
from which may be derived the mass function dAf / dm for the entire system of globulars in a given galaxy) in attempts 
to model this property of GCSs using our current understanding of present-day star formation in molecular clouds and 
their clumps. 

The work presented in this paper further explores the connections between the formation of stars, globular clusters, and 
galaxies, by concentrating on two of the primary attributes of GCSs: their specific frequencies (or total populations) and 
radial distributions. The next Section discusses the significance of these GCS properties as they relate to general questions 
on the formation of stellar clusters. In §3 below, data are compiled from the literature, and a geometrical deprojection 
algorithm applied, to construct volume density profiles for the stars, the X-ray gas, and the globular clusters in three 
galaxies with well studied GCSs: M87 (the cD at the center of the Virgo Cluster), M49 (giant elliptical in Virgo), and 
NGC 1399 (cD galaxy in the Fornax Cluster). An intercomparison of the density profiles of the three halo components 
yields a point-by-point measure of the GCS contribution to the total luminous mass in each system, as a function of 
three-dimensional gala-ctocentiic radius, Tgc. These GCS mass ratios, which are constant throughout each galaxy (at large 
Tgc) and the same in all three, are interpreted as an indication of the formation efficiency of globular clusters. In §4, these 
results are supplemented with broader considerations of the GCSs in other early-type galaxies and in the Milky Way halo, 
and with a direct estimate for the relative formation rate of open clusters in the Galactic disk, to argue for the existence 
of a universal efficiency, or probability, for the formation of a bound stellar cluster from a dense cloud of gas. Section 5 
discusses the implications of this for a few specific issues in galaxy and GCS formation and evolution, and §6 concludes 
with a summary. 

2. THE EFFICIENCY OF CLUSTER FORMATION FROM GLOBULAR CLUSTER SYSTEMS 

2.1. Cluster Form,ation Efficiencies 

The possibility that globular cluster systems can be connected not only to galaxy formation, but to ongoing star 
formation as well, is suggested by the fact that this latter process operates largely in a clustered mode: by mass, most new 
stars in the Milky Way appear in groups within the largest clumps in molecular clouds, and only rarely in true isolation. 
This is, by now, an established empirical fact (e.g., Lada et al. 1991; Zinnecker et al. 1993), and it is easily seen (Patel 
& Pudritz 1994) to be a direct consequence of the different power-law slopes in the mass function of molecular clumps 
{dAf /dm oc m~^-^, so that the largest clumps, which weigh in at 10^-10^ Mq, contain most of the star-forming gas mass in 
any molecular cloud) and the stellar initial mass function {dAf /dm oc m~^'^^, putting most of the mass in young stars into 
<^ 1 Mq objects). Similarly, young super star clusters, which tend to have sizes and masses that are generally comparable 
to the old globulars in the Milky Way, can account for as much as ~20% of the UV light in starburst galaxies (Meurer et 
al. 1995). 
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This is, however, not to say that all, or even most, stars are born into bona fide clusters that exist as coherent 
dynamical units for any length of time; on the contrary, it would appear that very few multiple-star systems actually 
form as gravitationally bound clusters. (A quantification of "very few" is one of the goals of this paper.) At some point 
during the collapse and fragmentation of a cluster-sized cloud of gas, the massive stars which it has formed will expel any 
remaining gas by the combined action of their stellar winds, photoionization, and supernova explosions. If the cumulative 
star formation efficiency of the cloud, SFE = M3tars/(M3tars + -M^gas), is below a critical threshold when the gas is lost, 
then the blow-out carries away sufficient energy that the stellar group remaining is unbound, and disperses into the field. 
The precise value of this threshold SFE depends on details of the dynamics and magnetic field in the gas cloud before its 
self-destruction, and on the timescale over which the massive stars dispel the gas; but various estimates place it in the 
range SFEdit ^ 0.2-0.5 (e.g.. Hills 1980; Mathieu 1983; Elmegreen & Clemens 1985; Verschueren 1990; Lada, Marguhs, 
& Dearborn 1984). 

A complete theory of star formation must therefore be able to anticipate the final cumulative SFE in any single piece 
of gas with (say) a given mass and density, and thereby predict whether or not it will form a bound cluster. No such 
theory yet exists. More within reach are considerations of a statistical nature, in which a distribution of cumulative SFEs 
(and, thus, the fraction of clouds that achieve SFE > SFEdit before self-destruction) is calculated for an ensemble of 
star-forming clouds, in either a protogalactic or a present-day context (see, e.g., Elmegreen 1983; Elmegreen & Clemens 
1985; Elmegreen & Efremov 1997). But even here, the systematic effects of a range in the initial masses of the gaseous 
clouds, or of changes in the assumed star formation law (which is probably the most outstanding "free" parameter here) 
have not been explored fully. It is therefore not clear that the dominant factor governing the formation of bound stellar 
clusters vs. unbound associations has been identified. 

It is nonetheless possible to constrain any forthcoming theories of star and cluster formation with empirical evaluations 
of the probability that a cluster-sized cloud of gas is able to achieve a cumulative SFE of at least the critical 20%-50%. 
This probability — which can equally well be interpreted as that fraction of an ensemble of massive star-forming clouds 
which manages to produce bound stellar systems — is referred to here as the efficiency of cluster formation. As was 
indicated above, an estimate of this will be taken first from a detailed investigation of the old globular cluster populations 
in three giant elliptical galaxies; general arguments will then be used to show that the result also applies to globular 
cluster formation in other galaxies, and to the formation of disk clusters today. 

A first (and very rough) assessment of the formation efficiency of globular clusters appeals to the most basic and 
most easily measured observable of a GCS: namely, its total population, Atot, summed over all cluster luminosities and 
galactocentric positions. This is usually compared to the integrated y-band luminosity of the parent galaxy's halo (as a 
measure of its unclustered field star population) through the so-called specific frequency (Harris & van den Bergh 1981): 

Sn = Mot X l0O-4(*^?+i5) oc Mot/iy,gai • (1) 

Even though any observational derivation of Sn necessarily involves (possibly large) extrapolations of data that directly 
count only the brightest clusters in some spatially limited subset of a GCS, these estimates tend usually to be surprisingly 
robust. To date, specific frequencies have been obtained for of order 100 galaxies, some 80% of which are early-type 
systems (see the compilations of Harris 1991; Durrell et al. 1996; Kissler-Patig 1997; Blakeslee 1997; Blakeslee, Tonry, & 
Metzger 1997; Harris et al. 1998). Given this bias in the observational database, some of what follows will make use of 
results that are really specific to elliptical galaxies; spirals are brought more explicitly into the discussion in §4 below. 

Specific frequency is easily related to the total masses of a galaxy's GCS and of its unclustered halo stars, if both 
the mean mass of the individual globulars in the system and the mass-to-light ratio of the overall stellar component are 
known. As is discussed further in §3.1 below, the basic distribution of globular cluster masses is quite similar from galaxy 
to galaxy, so that the mean cluster mass in the Milky Way GCS — (m)ci = 2.4 x 10^ Mq — serves as a reliable estimate for 
other systems as well. Meanwhile, the stellar mass-to-light ratios of large galaxies are essentially equal to the dynamical 
M/L in their cores^ which are baryon-dominated. As is also discussed below (§4.1), observations show that the core 
M/L of hot galaxies increases with luminosity; but, even though this turns out to be rather an important consideration 
in detail, it will suffice for now to use a single value of Ty = M/Ly = 7Mq Lq^ as a rough characterization of large 
ellipticals (van der Marel 1991; cf. Faber & Gallagher 1979; Binney & Tremaine 1987). The definition of Sn in equation 
(|^) may then be rewritten as (cf. Zepf & Ashman 1993; Harris & Pudritz 1994; Harris et al. 1998) 

S - 2500(^ ( ^ (2'\ 

^ \2Ax105MqJ [tMqLq^J Mst^,,- ^' 

To be clear, it is stressed again that the mass-to-light ratio applied here refers to measurements in the cores of ellipticals, 
and is not meant to include any dark matter other than stellar remnants. In any event, an inversion of this last expression 
yields 

- 4.0 X 10- 5„ ( .rJ^] ( " - 2 X 10- ( ^) , (3) 



M.tars ■ \2Ax105MqJ ^TMeLQ^ 

where Sn ~ 5 is suggested by the studies cited above as a representative value for early- type galaxies ranging in luminosity 
from My ~ —15 to My ~ —22. Of course, there is some uncertainty in this "mean" Sn, and some intrinsic scatter about 
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it. Moreover, the data suggest that Sn varies systematically with Ly^gai — the simple scaling A/tot oc Ly.gai that corresponds 
to a constant specific frequency is not strictly obeyed — particularly among dwarf ellipticals and between brightest cluster 
galaxies. These issues have been discussed in detail in the literature, and they are revisited here, in §2.2. With this caveat 
in mind, however, Mgcs/Mgtars ~ 2 x 10"'^ is useful as a crude, order-of-magnitude guide to a "typical" global GCS mass 
ratio in galaxy halos. 

The advantage of working in terms of this mass ratio is that, whatever its precise value in any one system at the current 
epoch, it is expected to be a fairly well conserved quantity that has not changed drastically over the course of GCS and 
galaxy evolution in a Hubble time. This is true even though the total population Atot of any GCS is bound to decrease 
over time, as individual globulars are disrupted by the cumulative effects of external gravitational shocks and internal 
two-body relaxation (see, e.g., Aguilar, Hut, & Ostriker 1988; Ostriker, Binney, & Saha 1989; Capuzzo-Dolcetta 1993; 
Murah & Weinberg 1997a, 1997b, 1997c; Gnedin & Ostriker 1997; Vesperini & Heggie 1997; but cf. Portegies Zwart et 
al. 1998). These dynamical processes discriminate preferentially against low- mass and low-density globular clusters, and it 
turns out that even if large numbers of such objects are lost, no great change is effected in the total mass of the remaining 
GCS. This conclusion rests on the fact that the mass functions of GCSs — particularly those in early- type galaxies — are 
shallow enough that most of the total Mgcs is contained in the most massive members of the system; but these are the 
least susceptible to dynamical destruction. The implication is that the total mass of a GCS may not change much over 
its lifetime; thus, if the current stellar mass of a galaxy were a reliable indication of the total amount of gas that was 
initially available to form stars, then the presently observed ratio Mgcs/Mstars ~ 2 x 10"'^ would serve as a rough, global 
estimate of the efficiency of globular cluster formation. Note that this reasoning also shows that galactic halos cannot 
be built up from disrupted globulars; see, e.g.. Ashman & Zepf (1992) and Harris & Pudritz (1994). (One of the main 
conclusions of this paper will be that, in fact, Mgtars is not always an adequate measure of an initial gas supply, and that 
Mgcs/Mgtars therefore cannot always be equated directly to a cluster formation efficiency. However, since most previous 
discussions of GCS specific frequencies do draw on this simple assumption at some level, it is of interest to develop some 
of its implications.) 

A simple numerical example should emphasize the robustness of Mgcs- Thus, consider a newly formed system of 
globulars whose masses are distributed according to a power-law mass function between some lower and upper limits: 

— K m~'^^ , me < m < mu , (4) 

init 

where 72 > and reasonable choices for m^ and m„ may be 500 Mq and 5 x 10^ M©. (The precise values are not critical 
here.) Now let this GCS be eroded over 10^" yr, keeping in mind that each cluster in the system has a definite, mass- 
dependent lifetime against the dynamical mechanisms mentioned above. In a rough approximation to a combination of 
complicated processes, there is then a mass scale (say m*) above which few or no globulars will have disappeared by the 
current epoch. Thus, the GCS mass function at m > m* is effectively unchanged, even after a Hubble time, from the 
initial distribution. However, clusters with initial masses m < m, can have been removed from the system, and the GCS 
mass function below may differ significantly from its initial form. If this low-mass end of the evolved d/V /dm can also 
be described by a power-law, then the full, present-day GCS mass function is just 

(dJV\ _ J K mj^ '^'^ m^^^ , me < m < 
dm)^^^ I Km~''^, m^ < m < mu , 

where 71 < 72 for consistency. 

Observed GCSs do, in fact, have mass functions that conform to roughly this sort of double power law (Harris & 
Pudritz 1994; McLaughlin 1994; McLaughlin & Pudritz 1996) — although there is somewhat more structure than this at 
the highest cluster masses in dAf /dm, it can be safely ignored for the purposes of the present argument. Empirically, 
the mass scale m* is always about 1.6 x 10^ Mq; at this mass, the logarithmic slope of dN /dm changes rather abruptly 
from 72 > 1 to 71 < 1, with precise values that can vary slightly from galaxy to galaxy There are strong arguments 
to support the notion that observed GCS mass functions at m > m* are indeed accurate reflections of the formation 
distributions (namely, the high-mass sides of dM /dm in observed GCSs show no significant variation with galactocentric 
radius inside a single system, and are quite similar — though not identical — from galaxy to galaxy; see Harris & Pudritz 
1994; McLaughlin & Pudritz 1996; Elmegreen & Efremov 1997). In addition, some theoretical studies (e.g., Surdin 1979; 
Okazaki & Tosa 1995; Vesperini 1997; Elmegreen & Efremov 1997) have suggested that the shallower slopes at lower 
masses might be caused entirely by dynamical evolution from an initial dM /dm oc to"'*'^, roughly as described above. In 
reality, this view could very well be too extreme: empirically, the system of young super star clusters in the Antennae, 
which is not likely to have been severely affected by dynamical destruction, may have a mass function of the form (|^) 
rather than (0) (Fritze-von Alvensleben 1999); and theoretically, such an initial shape to dM /dm may in fact be preserved 
during a Hubble time of dynamical evolution (Vesperini 1997). Thus, the intent here is not necessarily to argue for an 

■^It is this sharp change from an exponent of > 1 to one of < 1 at a well defined and nearly universal cluster mass that results in the 
characteristic peak of the globular cluster luminosity function in galaxies (the GCLF: the number of clusters per unit absolute magnitude, 
which is proportional to dV/rf In m; see McLaughlin 1994). This peak occurs at the luminosity corresponding to the mass mt, i.e., at 
My ~ —7.4, for the Tv — 2 Mq typical of globulars; its universality is what makes the GCLF attractive as a distance indicator (e.g., 
Jacoby et al. 1992; Whitmore 1997). 
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explanation of present-day GCS mass spectra in purely evolutionary terms, but to show that even if such is adopted as 
an extreme possibility, the implied time dependence of total GCS masses is quite weak. 

Given equations (H) and (g), then, straightforward integrations allow for comparisons of the present GCS population 
to the initial one: 



"tot 



J\f°^,- _ fmuV' ' ri - (m*/TO„)^^-i + [(72 - 1)/(1 -7i)][l - (m,/m*)i-^i] 



,m*y [ {mu/mep^-'^ -1 

and of the current total mass to its initial value: 

M°^i fmA^-''' r(m„/TO*)2-^^_l + [(2-72)/(2-7i)][l-("^^M*)'"^^: 



M^cs \muj I 1 - (m£/m„)2-72 



(6) 



(7) 



Values of 71 — 0.4 and 72 = 1.7 may be taken as representative of the typical situation in elliptical galaxies (Harris 
& Pudritz 1994; McLaughlin 1994, 1995; McLaughlin & Pudritz 1996). With me = 500 Mq, m, = 1.6 x lO^M©, and 
m„ = 5 X 10^ Mq, it is then easily seen that 

0.036 and ^£^=0.76 (8) 

in this simplistic scenario. Thus, even though dynamical evolution reduces the total population of this hypothetical GCS by 
more than an order of magnitude, it diminishes the total mass by only ~25%. (The same line of argument also shows that 
the often appreciable uncertainties in the exact form of dJ\f /dm at cluster masses m < — or even blind extrapolations 
from observations of only brighter clusters — do not seriously bias observational estimates of Mgcs in any galaxy.) 

This result is best interpreted as an estimate for the possible depletion of a GCS globally (i.e., spatially averaged over the 
whole of a parent galaxy), because the values of 71, 72, and m* used are taken from observations of dAf /dm in projection, 
often over fairly large areas on the sky. Thus, the normalization K in equations (jj) and (||) can be taken to include an 
implicit integration over galactocentric radius, and the exponents 7 and mass scale to* to be appropriate averages. The 
importance of these considerations is not actually clear, since radial dependences are not expected theoretically in current 
models for the initial GCS mass function (McLaughlin & Pudritz 1996; Elmegreen & Efremov 1997) and observational 
evidence for such trends at the current epoch is slight to non-existent (e.g., McLaughlin et al. 1994 and Harris et al. 1998; 
though see also Kavelaars & Hanes 1997 and Murali & Weinberg 1997b). Nevertheless, it should be noted that dynamical 
friction, gravitational shocks, and even evaporation couple to galactic tidal fields in such a way that if a GCS's present-day 
mass function were sculpted to any great extent by these processes, then the local ratios M°^^ / Nl'^^ and M °^g^/Mg"g* would 
be expected to increase outwards from the center of the system, from values less than equation (ph at small galactocentric 
radii to near unity at large distances (e.g., beyond a stellar effective radius; see especially Murali & Weinberg 1997b). 

Returning now to the main issue, the global efficiency of globular cluster formation in a galaxy may be defined as 

6ei ^ M^^/M-^: , (9) 

where M™g* refers to the total gas supply that was available to the protogalaxy, whether in a rapid monolithic collapse or, 
more likely, in a slower assembly (at high redshift) of many distinct, subgalactic clouds. From the preceding discussion, it 
is clear that if a newly formed GCS were able to evolve "passively" — that is, if its mass were potentially vulnerable only 
to reductions caused by forces internal to an isolated parent galaxy — then a present-day observation of Mgcs would suffice, 
to at least a ~25% level of accuracy globally, as a measure of Mg"*. Passive evolution also implies mass conservation 
for the galaxy as a whole, so that a reliable value for Mg"^* could be obtained, at any epoch, from an inventory of the 
total mass in gas plus unclustered stars and stellar remnants. Thus, a good observational estimate of the global cluster 
formation efficiency in this simple situation is provided by 

-(Wgas + JWstars 

which is expected to be nearly a time-independent quantity]^ 

Realistically, of course, galaxies are not such pristine entities; interactions between systems may effect significant changes 
in both Mgcs and (Mgas + Mstars) over time, making it difficult to evaluate (or even clearly define) initial GCS and gas 
masses. For example, the formation history of an elliptical may include major mergers with gas-rich spiral systems, and 
it has been proposed that these might involve the formation of new globular clusters and substantial increases in Mgcs 
(e.g.. Ashman & Zepf 1992; Kumai, Basu, & Fujimoto 1993a). Alternatively, additional globulars might be captured by a 
large galaxy that accretes small, gas-poor systems (e.g.. Cote et al. 1998) or accumulates tidal debris from other members 

^To be fully consistent, of course, the observational estimate of cluster formation efficiency should read e^i = Mgcs/(Afgas 4- MstniB + -Afgcs); 
but, as §§3 and 4 below will show, Mgce turns out always to be so small, relative to (A/gas + Mstars), that its contribution to the denominator 
of Eel is utterly negligible. 
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in a cluster of galaxies (see Muzzio 1987). However, processes such as these also increase the total masses of gas and/or 
stars in the final system. In large galaxies (i.e., in those that do not suffer significant gas loss from strong galactic winds; 
cf. §4.2), any changes in the ratio Mgcs/(Mgas + Mstars) are bound to be smaller than changes in either of these quantities 
individually. The presently observed eci in any large galaxy with a complex dynamical history is then a mass-weighted 
average of the true ed in some number of discrete star-formation events and individual accreted systems; and in the limit 
that the underlying efficiency of cluster formation is independent of environment and epoch (as the results of §§3 and 4 
will suggest is nearly correct), it is still true that eci — Cci — ^just as in the case of passive evolution, and subject only to 
the small error associated there with the effects of CCS erosion by cluster evaporation and tidal disruption. 

In gas-poor elliptical galaxies and the halos of spirals, equation ( [lO| ) reduces to eci — -Mgcs/Mstars, so that — as was 
mentioned above — the cluster formation efhciency can be obtained directly from obervations of the global CCS specific 
frequency (cf. eq. More generally, however, ellipticals can contain large amounts of hot. X-ray emitting gas (e.g., 

Forman, Jones, & Tucker 1985), and when it happens that Mgas is non-negligible relative to Mgtars, measurements of 
Sn are only low-order approximations to EcI (which is itself only an estimate of the real efficiency ed). Keeping this in 
mind, the observed systematics of both global and local CCS specific frequencies provide a useful starting point for an 
examination of the extent to which cluster formation efficiencies might vary from galaxy to galaxy, or from place to place 
within a single system. 

2.2. The Specific Frequency "Problems" 

As was intimated in the discussion around equation (|^), not all galaxies have the same globular cluster specific frequency. 
There are several aspects to this: (1) Spiral galaxies tend to have lower Sn — at least when their CCS populations are 
normalized to their total (bulge plus disk) light — than do average ellipticals (Harris 1991). (2) Although a constant 
Sn « 5 is a fair approximation for many normal giant ellipticals, fits to the data suggest that, in fact, Sn oc iygai"'^ 
(Kissler-Patig 1997; cf. Santiago & Djorgovski 1993). At the extremes of the early-type galaxy sequence, moreover, (3) the 
Sn values for centrally dominant or first-ranked members of galaxy clusters show a much steeper dependence on galaxy 
luminosity — Sn csc gai or so is implied by the correlations of Blakeslee (1997) and Harris et al. (1998) — and (4) specific 

frequencies actually appear to increase towards lower luminosities in dwarf ellipticals: Sn oc -^ygai (Miller et al. 1998; 
see also Durrell et al. 1996). Thus, some brightest cluster galaxies, and some of the faintest dwarfs, have global specific 
frequencies that exceed the "average" by factors of 3 or 4. All in all, observed Sn values span a range from ~1 to ~20n 

As will be shown in §4 below, the first two of these points can be understood fairly easily. Here, it suffices to note that (1) 
when only the halo stars and globular clusters in the Milky Way are considered, and allowances made for the lower Pop. II 
mass-to-light ratio in our Galaxy vs. large ellipticals, Cci here is no different from that in "normal" early-type systems. 
It is reasonable to expect that this result — which is contrary to the claim made by Zepf & Ashman (1993) — should hold 
in other spirals as well. Also, (2) the observed dependence of Sn on Lv,gai among non-brightest cluster ellipticals is 
fully consistent with that expected from equation (0) for gas-poor systems with a uniform efficiency of cluster formation 
(cci — Mgcs/Af stars) and a mass-to-light ratio (Ty) that scales with luminosity in the way observed for galaxies in the 
fundamental plane. However, if taken at face value, items (3) and (4) just above appear to imply that GCS formation 
efficiencies can vary by as much as an order of magnitude from galaxy to galaxy, in a non-monotonic fashion and for 
reasons that have never been understood. These issues are not so readily dismissed. The high Sn "phenomenon" in 
many first-ranked cluster ellipticals, especially, has been discussed at length in the literature (e.g., McLaughlin et al. 1994; 
Harris et al. 1995; West et al. 1995; Forbes et al. 1997; Blakeslee 1997; Blakeslee et al. 1997; Harris et al. 1998), and is 
referred to here as the first specific frequency problem. Equally non-trivial, though only more recently recognized, is the 
inverse variation of Sn with Ly^ged among dwarf ellipticals. 

Local specific frequencies, which are obtained for a single GCS by comparing its projected radial distribution to its 
galaxy's surface brightness profile (normalized as in eq. Q), are also of interest here; they are related in the obvious way 
to the efficiency of cluster formation as a function of (projected) galactocentric radius, Rgc- Of particular importance 
is the fact that Sn is seen to be an increasing function of Rgc in many galaxies. This is equivalent to the well known 
fact that GCSs are often less centrally concentrated than the stellar halos of their parent galaxies. This statement has 
two contexts, which must be separated. First, it seems that the GCS of most any galaxy — spiral or elliptical, small or 
large — has an apparent core radius that is significantly larger than that of the unclustered field stars (e.g., Djorgovski 
& Meylan 1994; Lauer & Kormendy 1986; McLaughlin 1995; Forbes et al. 1996; Minniti, Meylan, & Kissler-Patig 1996; 
Cote et al. 1998); that is, the central regions of GCSs show projected radial distributions with slopes that are significantly 
shallower than those of the light profiles on kpc scales, so that the density ratio of the two halo components is an increasing 
function of radius. Second, the GCSs of some large galaxies are spatially more extended than the stellar halos on larger 
scales, Rgc ~ 10-100 kpc (Harris & Racine 1979; Harris 1986). However, this is a less general result than that relating to 
the large core radii of globular cluster systems. It is true of many dominant cluster galaxies like M87 in Virgo (see the 
discussion in McLaughlin et al. 1995), and also of the M31 halo (Pritchet & van den Bergh 1994); but many elliptical 
systems (e.g., Durrell et al. 1996; Kissler-Patig et al. 1997), and the Milky Way itself (Harris 1976; Zinn 1985), show 
GCS and halo-star distributions that are in quite close agreement (and, thus, a constant local Sn) outside of any core 
regions. The physical factors that govern whether the GCS density falls off at the same rate as or more slowly than the 

has also been claimed (Harris 1991; West 1993; Kumai, Hashi, & Fujimoto 1993b) that early-type galaxies in sparse groups have 
systematically lower Sn than those in high-density cluster environments. However, this correlation is uncertain, and apparently rather weak; 
if it does prove to be important, it will most likely be as a higher-order correction to the more basic trends discussed here. 
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unclustered halo density in a given system have not yet been identified. (No case has ever been found in which a GCS 
radial distribution is steeper than a galaxy's light profile.) 

The first of these points is not easily related to the radial dependence of cluster formation efficiencies because, as was 
suggested in §2.1, the dynamical mechanisms that act to destroy globulars work most efficiently at small Rgc, where the 
stellar densities are highest. Although the cores of GCSs could quite plausibly have always been somewhat larger than 
those of halos or bulges generally (e.g., Harris & Pudritz 1994; Cote et al. 1998), it seems inevitable that any initial 
discrepancies must have been augmented by an inside-out erosion of the cluster system from tidal shocking, dynamical 
friction, and evaporation (Ostriker et al. 1989; Capuzzo-Dolcetta 1993; Murali & Weinberg 1997b). Detailed modelling is 
therefore required to disentangle the effects of formation and evolutionary processes in this instance. At galactocentric 
radii greater than a few kpc, however, the timescales on which the evolutionary mechanisms operate can exceed a Hubble 
time, and the studies cited here show that dynamics will not significantly affect GCS radial distributions beyond about one 
effective radius in the stellar halo (see also Aguilar et al. 1988). Thus, to the extent that the assumption eci — Afgcs/Mgtars 
is valid, the distension of GCSs relative to stellar halos and the concomitant systematic increase in local Sn in the outer 
parts of some galaxies would seem to say that bound clusters were sometimes (but not always, and for unknown reasons) 
more likely to form at larger i?gc, in gas that was presumably at lower ambient densities and pressures. This situation 
will be referred to as the second specific frequency problem. 

It is clearly of interest, from the points of view of both galaxy and star formation, to verify whether these global and 
local Sn variations really do reflect differences in Cci — and, if so, to understand their origins. Accordingly, much attention 
has been paid to the first Sn problem in particular, usually in efforts to empirically correlate the Sn measurements for 
brightest cluster galaxies (BCGs) with global properties of the galaxies themselves or of the clusters they inhabit. One 
such correlation does appear to be of fundamental importance, as will be discussed momentarily, but it suggests that Sn 
varies in the face of constant Cd- Indeed, no proposed solution of the first specific frequency problem has identified a 
physical mechanism that could possibly lead to the superefficient production of bound star clusters in some galaxies (see 
the references given above). In fact, one suggestion — that high global specific frequencies (meaning Sn ^ 5) in centrally 
located BCGs result from the acquisition of globulars and stars which are tidally stripped from the outer, high- at regions 
of other galaxies in a cluster — actually appears to rely on the existence of the second specific frequency problem (see 
Forbes et al. 1997; Cote et al. 1998); but this is itself essentially unexplained.^ Nor, to date, has any explanation (or 
prediction) been given for the high Sn that have now been observed in many faint dwarf galaxies. 

The view is taken here that the second specific frequency problem is inseparable from the first, in the specific sense that 
both relate to the possible influence of environment on the probability of bound cluster formation. Any real insight into 
either one of these issues must therefore include a better understanding of the other. This view is in keeping with recent 
evidence (Forbes et al. 1997; Kissler-Patig 1997; van den Bergh 1998) for at least a rough correlation between global Sn 
values and GCS density-profile slopes in elliptical galaxies: those GCSs with the highest specific frequencies tend also to 
have the shallowest radial distributions. Such a correspondence hints strongly at (although it does not directly show) an 
intimate connection between an increasing local Sn and a high global Sn- 

A significant development in this field has come with the recent recognition of a correlation between the global specific 
frequencies of BCGs and the X-ray luminosity of the hot gas in their parent clusters (West et al. 1995; Blakeslee 1997; 
Blakeslee et al. 1997). A relatively large and homogeneous sample of BCGs shows Sn values that increase from about 4, 
in clusters with X-ray luminosities (on 500-kpc scales) Lx ~ 10'*^ erg s""'^ at 0.5-4.5 keV, to Sn — 12 for Lx ^ erg 
s~^. Three different interpretations have already been attached to this result: 

West et al. (1995) propose that it points to the existence of intergalactic globulars which collect naturally at the centers 
of galaxy clusters and are superimposed in projection on any galaxy there. This would certainly boost the apparent total 
population A/tot of a central BCG, so that globulars need not have formed in situ with abnormally high specific frequency; 
but they would still have had to do so elsewhere in the cluster (see also Cote et al. 1998; Harris et al. 1998). 

In a different approach, Blakeslee (1997) and Blakeslee et al. (1997) consider Lx to be mainly a reflection of the 
temperature of intracluster gas, and thus a measure of a cluster's mass and the depth of its potential. They then 
claim that the observed Sn-Lx correlation is consistent with the globulars in BCGs having formed in numbers directly 
proportional to the total mass, including dark matter, of their entire galaxy clusters (see Kavelaars 1998 for a similar 
suggestion on the scale of individual galaxies). In this scenario, cluster formation would obey a universal efficiency of 
sorts — although one which is quite different from that discussed here — and the first Sn problem would result from some 
process that caused the BCGs in more massive clusters to be systematically underluminous (i.e., underefficient in the 
production of unclustered field stars). Here, however, the role that the dark matter content of a galaxy cluster might play 
in the formation of globular clusters is not clear. Moreover, any rule such as Af™*/(M™g* -I- Mdark) — constant clearly 
can apply over only a limited range of conditions. 

^It should be noted that the second ^jv problem, and local specific frequencies generally, are not typically discussed in terms of cluster 
formation efficiencies, even though it is well understood that global Spf are of interest in this regard. The reason for this may be that mismatches 
between GCS and stellar radial distributions in E galaxies are usually attributed either to the globulars having formed "slightly before" the bulk 
of the unclustered halo stars (in a dissipational collapse scenario; Harris 1986), or to the expectation that a cluster system will be "puffed up" 
relative to the final light profile of an elliptical that forms in one or more major mergers of gas-rich systems with pre-existing GCSs (Ashman & 
Zepf 1992; Zepf & Ashman 1993). In either case, no particular radial variation in e^i would necessarily be required. However, the first of these 
arguments has not been sufficiently well developed to account quantitatively for the observed magnitude of the second Spf problem; and the 
second is actually inconsistent with the observation that systems with higher global Sm have spatially more extended GCSs (see the discussion 
of Forbes et al. 1997). Nor is it apparent, from either these suggestions, why there are also many galaxies whose GCSs and stellar halos have 
the same density profiles. 
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Thus, Harris et al. (1998) instead pursue the hypothesis that GCSs formed with a constant defined essentiahy as 
in equation above. This is adopted as a working assumption by Harris et al., who construct plausibiHty arguments 
for a universal value of ed ~ 10~^ (as could be expected) and suggest that the higher global Sn in more X-ray luminous 
BCGs might have resulted from stronger (partial) galactic winds which heated larger fractions of their initial gas supplies 
to keV temperatures after an early bout of star and globular cluster formation. Such hot gas would be largely removed 
from further star formation — but still confined to the vicinity of these central cluster galaxies — and the final BCGs would 
be optically under luminous as a result. 

It is this basic idea that is closest to the discussion in this paper (although many aspects of the quantitative analysis 
here differ significantly from those in Harris et al. 1998). That is, the fundamental point of the global Sn~Lx correlation 
for BCGs is the associated implication that the first specific frequency problem may stem, not from intrinsic variations 
in the global ed of these systems, but from the fact that they can contain large amounts of gas. High gas mass ratios 
invalidate the assumption that Cci — -A^gcs/Afstars, and therefore cause Sn to overestimate the true globular cluster 
formation efficiency. If it is correct to infer from this that the "excess" globulars in systems with high global Sn are 
associated with hot gas, then the fact that the distribution of gas is generally more extended than that of the stars in 
ellipticals (e.g., Trinchieri, Fabbiano, & Canizares 1986; Fabbiano 1989) suggests a possible explanation for the second 
specific frequency problem as well as the first. These two issues would then be closely related — possibly even just different 
aspects of a single phenomenon — as was suggested above; and both might be settled simply by measuring the full estimator 
ed = Mgcs/{Mgas + -Mstars) both locally and globally in individual systems. 

In the next Section, this universal-ed hypothesis is tested, for the first time, by evaluating the observable Cd as a function 
of galactocentric radius in each of the galaxies M87, M49, and NGC 1399. It turns out that the inclusion of the X-ray 
gas in M87 is indeed key to the alleviation of both the first and the second Sn problems there: once this is done, Cci is 
constant locally, and essentially the same globally, in all three of these systems. The discussion is expanded to include 
other giant ellipticals and BCGs, as well as dwarf ellipticals, in §4. It is shown directly that all available CCS data are 
consistent with a single efficiency of cluster formation. Notably, the increase of Sn towards low Ly.gai in the dwarfs (Miller 
et al. 1998; Durrell et al. 1996) can be explained in terms of this universal Cd once allowances are made for the effects of 
the violent gas outflows (a more extreme version of the feedback envisioned by Harris et al. 1998 for BCGs) that were 
driven by early bursts of star formation in these low-mass systems. 

3. CASE STUDIES OF THREE GALAXIES 

M87, M49, and NGC 1399 are ideally suited to a closer examination of these issues, for three reasons: (1) They are 
relatively nearby, making it possible to observe internal trends within the systems in the first place. M87 and M49 both 
are assumed here to lie at a distance oi D = 15 Mpc, appropriate for the core of the Virgo Cluster (Pierce et al. 1994; 
Freedman et al. 1994). A relative distance modulus of A{m — A/) [Fornax — Virgo] — 0.2 (e.g., McMillan, CiarduUo, & 
Jacoby 1993; Tonry et al. 1997) then puts NGC 1399 at D ~ 16.5 Mpc. (2) Each is a giant elliptical (M87 and NGC 1399 
are actually type cD), and is therefore accompanied by a large GCS: A/tot — 13 500 for M87; 6 900 for M49; and 4 700 
in NGC 1399 (e.g., Harris 1991; see also §5 below). And, most importantly, (3) even this small sample of three galaxies 
suffices to bring out all aspects of the two specific frequency problems summarized in §2.2: M87, as is well known, suffers 
from both; M49 does not show the first, but has long been thought to exhibit the second; and NGC 1399 has been claimed 
as an example of the first problem, but not the second. 

Because of their proximity, these systems have been well studied in all respects of importance here. _B-band surface 
photometry extends in M87 to a projected galactocentric radius of about 22' ~ 95 kpc (de Vaucouleurs & Nieto 1978, 1979); 
in M49, to a similar distance (King 1978; Caon, Capaccioh, & D'Onofrio 1994); and in NGC 1399, to i?gc ~ 14f5 = 70 
kpc (Killeen & Bicknell 1988; Caon et al. 1994). These measurements are easily converted to projected stellar mass 
density profiles Estais(^gc), as mass-to-light ratios have been obtained from a spectroscopic analysis applied uniformly to 
all three systems (van der Marel 1991). X-ray observations, which have been made with several different satellites, yield 
gas densities over still larger areas around each galaxy. Here, use will be made of the RO SAT data for M87 (taken from 
Nulsen & Bohringer 1995); of i?OS'^T data for M49 (Irwin & Sarazin 1996); and of Einstein observations in NGC 1399 
(Killeen & Bicknell 1988). These studies all give the volume electron density ne(rgc), as a function of three-dimensional 
galactocentric radius, derived either from fits of a parametric function to the projected X-ray brightness profile or from a 
direct deprojection of the same. Finally, integrated optical magnitudes have been measured for globular clusters in these 
galaxies at radii Rgc ^ 1 — 100 kpc — generally, over most of the areas covered by the surface photometry of unresolved 
halo stars. These are used to construct projected GCS number densities, Nc\{Rgc), which are trivially converted to mass 
profiles, Ed, once a mean globular cluster mass is specified. For M87, relevant investigations are those of Harris & Smith 
(1976), Harris (1986), Grillmair, Pritchet, & van den Bergh (1986), Lauer & Kormendy (1986), Cohen (1988), McLaughlin 
et al. (1993, 1994), McLaughhn (1995), Whitmore et al. (1995), and Harris et al. (1998). Data on M49 can be found 
in Harris & Petrie (1978), Harris & van den Bergh (1981), Harris (1986), Harris et al. (1991), and Lee, Kim, & Geisler 
(1998). The GCS of NGC 1399 has been observed by Hanes & Harris (1986), Bridges, Hanes, & Harris (1991), Wagner 
et al. (1991), and Kissler-Patig et al. (1997). 

While surface brightness profiles and X-ray gas densities often can be obtained for most of an entire galaxy with a 
single set of observations — as is the case in the studies cited above — this is generally not so for GCS radial distributions. 
Older, photographic surveys, which did have large {R ~ 100 kpc) fields of view, were hampered by rather bright limiting 
magnitudes that became systematically more so towards smaller galactocentric radii. Newer, CCD photometry achieves 
deeper limiting magnitudes, and digital data reduction techniques can quantify and correct for any spatial variations; but 
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until very recently, CCD fields of view were quite small, and best suited to analyses of the central regions {Rgc < 10 kpc) 
of GCSs. Thus, the first step in comparing the stellar, gaseous, and globular cluster contents of any large galaxy is to 
combine existing photographic and CCD CCS data into a single, composite radial distribution that gives the projected 
density -/Vci(^gc)j to a uniform limiting magnitude, over as large a range of galactocentric radius as possible. This is done 
for each of M87, M49, and NGC 1399 in §3.1. But then, the surface densities of stars and GCSs are related to their 
volume densities by an integral of the form E oc / pdz, while the observed X-ray brightness along some line of sight is 
proportional to / p^^^ dz. A direct intercomparison of the raw GCS, optical, and X-ray imaging data in a galaxy therefore 
makes little sense. Thus, to discuss these three components simultaneously in M87, M49, and NGC 1399, either (1) their 
GCS and stellar surface densities can be deprojected and compared to their derived gas volume densities, or (2) the gas 
surface densities, Sgas(-Rgc) « J Pgasirgc) dz, can be obtained from the published volume density profiles and compared to 
the directly observed Sstars(^gc) and Sci(-Rgc). The results of these two separate exercises are presented in §§3.2 and 3.3. 

3.1. GCS Surface Densities 

When joining GCS radial distributions from two or more separate studies of the same galaxy, two important factors 
must be taken into account. First, the GCS of a distant galaxy (essentially, any with D ^ 5 Mpc) is usually identified 
as a centrally concentrated excess of point-like sources, above a more uniform distribution of unresolved background 
galaxies and foreground stars. The surface density of this "background," Nt,, is subtracted from the total density of point 
sources at any Rgc to form the surface density profile Nc\ of the GCS alone; but Nj, is a function of seeing conditions, 
detector characteristics, and limiting magnitude. Background corrections must therefore be applied separately to each of 
the datasets to be combined. Second, differences in limiting magnitude (call this Vnm) are of further importance because 
a deep survey will clearly see more globulars, at a given -Rgc, than one with a bright limit. Before they can be simply 
matched onto one another, then, distinct background-corrected profiles Nci{Rgc) must be normalized to the same effective 

Vlim- 

The published GCS observations that will be used here for M87, M49, and NGC 1399 are summarized in Table |l|. These 
papers were selected from the more numerous references cited above, simply because they tabulate all of the relevant data 
in full. The third column of Table gives the minimum and maximum Rgc contained in each study's field of view (note 
that 1' = 4.363 kpc for a distance of D = 15 Mpc to M87 and M49, while D = 16.5 Mpc for NGC 1399 imphes 1' = 4.800 
kpc); Column 4 gives the background surface density to the limiting magnitude in Column 5; and Column 6 lists the 
scalings applied to each dataset to bring them all to a common Vnm- 

Generally, the background levels have been determined and applied by the original authors, either from observations 
of a blank field near the observed GCS, or from "on-frame" extrapolations of counts at the edges of the GCS field itself (see 
Harris 1986); these values are reproduced in Table 0only for completeness. The single exception to this is the M49 study 
of Lee et al. (1998), for which Nt has been derived here, using the method of Harris (1986). (It is possible, particularly 



Table 1 

Published GCS Radial Distributions 



Galaxy and 


Source 


Radial Coverage 


Background Density 




Scale 


GCLF Parameters 




(arcmin) 


(arcmin"^) 


(mag) 


Factor 


M87 


McLaughlin (1995)^ 


0.153 - 1.965 


3.5 ±0.8 


23.9 


1.000 


V° = 23.70 cry = 1.39 


McLaughlin et al. (1993) 


1.210-9.090 


6.3 ±0.4 


24.0 


0.954 




Harris (1986) 


1.000- 22.52 


5.8 ±0.3 


23.4 


1.35 


M49 


Harris & Petrie (1978)^ 


2.000- 19.69 


4.3 ±0.3 


23.3 


1.000 


V" = 23.75 cry = 1.30 


Harris & van den Bergh (1981)^ 


0.732 - 12.15 


3.3 ±0.5 




0.86^ 




Harris et al. (1991) 


0.600- 2.304 


2.88 ±0.45 


23.8 


0.707 




Lee et al. (1998) 


0.900- 7.167 


2.95 ±0.30 


23.45 


0.892 


NGC 1399 


Kissler-Patig et al. (1997) 


0.761 -9.463 


6.1 ±0.3 


24.0 


1.000 


V° = 23.90 cry = 1.20 


Hanes & Harris (1986) 


0.734- 13.97 


5.40 ±0.45 


23.4 


1.575 



^See also Harris et al. (1998). 

^As given by Harris (1986). 

■^See the discussion of Harris (1986). 
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with on-frame estimates, to overestimate Nh and produce spurious curvature at large radii in a subtracted A^ci profile. As 
is discussed below, however, this is not a major concern in the composite GCS radial distributions obtained here.) The 
limiting magnitudes Mim in Table |l] are also those quoted in the original papers.]^ These have been used to renormalize the 
individual A^ci profiles, by adopting — for computational convenience — a standard Gaussian approximation to the globular 
cluster luminosity function (GCLF, dAf/d In to; cf. §2.1). That is, with the number of globulars at a given apparent 
magnitude given roughly by dAf /dV oc exp[— (y — y°)^/2cry] (see Harris 1991), then the density of clusters to some 
Hmiting magnitude, say Nci{V < Viim,2), is referred to the density up to some other Umit, i.e., Nci{V < Mim,i); through 

iVci(V^ < Mim.i) ^ 1 + erf [(Mima ^'^'°)/V2ay] 
Na{V<V,i^,2) 1 + erf [(l^i„,^2-^°)/V2ay] ' 

The right-hand side of this expression is just the scaling factor given in Column 6 of Table |l|, where Vlim,! for each 
GCS has been taken as the observed limiting magnitude in the first of each set of references. Thus, the published radial 
distributions for M87 have all been normalized to an effective T4im = 23.9; for M49, to T4im = 23.3; and to Vlim = 24.0 
in NGC 1399. The Gaussian peaks V° and dispersions ay used in equation (|l^) are given in the first column of Table 
1^. These parameters are taken from the fits of Whitmore et al. (1995) and Harris et al. (1998) to the M87 GCLF, from 
Lee et al. (1998) for M49, and from Kohle et al. (1996) for NGC 1399. Note that all of the scale factors derived here are 
rather modest, and thus do not introduce significant uncertainty into the final, combined Nd profiles. 

A few practical notes should be made regarding this scaling procedure. First, it is applied to the GCS densities at 
every observed galactocentric radius, and it therefore assumes that there is no intrinsic radial variation in the GCLF or 
the underlying GCS mass function. As was also mentioned in §2.1 above, this is a reasonable assumption that is directly 
supported by the observations themselves. Second, the limiting magnitudes for the photographic studies of Harris & 
Petrie (1978) and Harris (1986) are rather poorly determined; Harris quotes only approximate ranges for Viim- The values 
adopted here, which are fully consistent with these ranges, have been chosen to give scaled Nc\{Rgc) that agree well with 
more recent CCD studies with accurately known Vim- And third, the scaling applied to the M49 GCS counts of Harris 
& van den Bergh (1981) has not been calculated in this way at all; rather, the factor of 0.86 is derived, by Harris (1986), 
by requiring a good overall agreement with the A'ci profile of Harris & Petrie (1978). 

The results of this exercise are illustrated in Fig. |[ Clearly, the different projected density profiles in each GCS 

®Not all of the photometry in these studies is in the V band, however. A mean globular cluster color of V — R = 0.5 has been adopted to 
convert the R-band data of McLaughlin (1995) to V magnitudes; a mean V — Ti = 0.45 brings the Washington photometry of Lee et al. (1998) 
onto the same system; and approximate transformations between photographic and V-band magnitudes are discussed in Harris' earlier papers 
on M87 and M49. 
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corroborate each other. It remains only to extract a single A^ci at each Rgc, in each system, from the overlapping data 
that have been collected. It is not generally possible to simply average different densities at a given Rg^, since profiles 
such as these are constructed by counting clusters in a series of wide, circular annuli that need not coincide from study 
to study. Instead, Fig. |l| shows that a perfectly acceptable approach is to simply adopt one of the (scaled) profiles over 
some range in i?gc, and then switch to any other that might extend beyond that range. At this point, it must be kept 
in mind that ultimately the data will be deprojected to yield volume GCS densities, nci{'^gc)- The method developed in 
the Appendix to do this requires that the inner and outer radii be known for the annuli used to define the composite Nd- 
This further restricts the useful published datasets to those which contain this information on the radial bins as well (that 
is, in addition to just a density and a mean radius for each annulus). 

In M87, then, a composite GCS profile is derived by using the data of McLaughlin (1995) from Rgc — 0.'153 to 
Rgc = 1.'870; those of McLaughHn et al. (1993) from 1.'870 to 7f870; and those of Harris (1986) from 7.'870 to 22'.52. 
(In order to make annuli that do not overlap in the final iVd, some interpolation is applied to the datasets around 
the matching points at 1.'870 and 7.'870.) For M49, the scaled profile of Harris & van den Bergh (1981) is used over 
0f732 < Rgc < 10^96, and that of Harris & Petrie (1978) for 10^96 < Rgc < 19f69. And in NGC 1399, the densities of 
Kissler-Patig et al. (1997) are adopted for 0f761 < Rgc < 7^571, with the data of Hanes and Harris (1986) continuing 
on out to Rgc = 13f97. In all cases, CCD data have usually been chosen over photographic data in any regions of 
overlap. However, superseding this preference is a concern over possible errors in the background estimation Nt; where 
this seems too high — for example, when the background-subtracted Nc\ at the largest radii in a smaller-scale CCD study 
drop significantly below the densities indicated by a more extensive photographic survey (as in the outermost data points 
of Kissler-Patig et al. 1997 for NGC 1399) — the photographic profiles are given precedence. Enforcing consistency in this 
way between independent GCS datasets ensures that the question of potentially overestimated Nh values remains an issue 
only at the largest radii (beyond the regime of any overlap) in the widest-field observations of each galaxy. 

The composite radial distributions thus obtained are listed in the third columns of Tables ^ through ^, as average 
surface densities in series of concentric circular annuli (with inner, outer, and average projected radii given in the first 
and second columns of the Tables). They are also drawn as the bold lines in the three panels of Fig. H^. There are two 
details of particular note. First, the densities for the innermost two annuli in M49 have been scaled up from the (already 
renormalized) GCS profile of Harris & van den Bergh (1981), by a further factor of 1.1, to put the composite A^ci at small 
radii into better agreement with the CCD data of Harris et al. (1991) and Lee et al. (1998). This ad hoc correction most 
likely just reflects some residual incompleteness in the photographic GCS counts at small Rgc, and in any event does not 
significantly affect the results which follow. Second, all of the Nc\ profiles have been smoothed somewhat by rebinning, 
in order to allow for a more stable deprojection in §3.2 below: the outermost six annuli of Harris (1986) are combined 
into one for the M87 profile; the outermost 11 rings of Harris & Petrie (1978) are blended into 4 for M49; and the last 
8 radial bins of Hanes & Harris (1986) become 4 in the NGC 1399 GCS. (As an added benefit of this smoothing, any 
possible errors in the subtracted Nj, are expected to affect only the outermost one or two points in the final Nd profiles.) 
Similarly, counts in the "Northeast" and "Northwest" sectors of NGC 1399 in the analysis of Kissler-Patig et al. (1997; 
see their Table 8) have been added together in every annulus they define; and these annuli have been combined in pairs 
to give the coarser profile in Table ^ here. 

Given these radial distributions, it is a simple matter to derive the total mass surface density profiles of the GCSs. 
This requires, first, another multiplicative scaling of 7Vci(i?gc) in each case, so that the number density of all globulars is 
properly reflected. (Note that the observations used here reach limiting magnitudes which are just at, or brighter than, 
the peak of the symmetric Gaussian that best fits the GCLF. This implies that most of the faint globulars are lost in the 
noise in these studies, and that half or fewer have been directly observed at any Rgc-) These final correction factors are 
obtained from equation ( |Tl] ) above, by setting Viim,i = oo and putting Viun,2 at the effective limit for each composite A^ci- 
In this way, it is found that the M87 distribution in Table || should be multiplied by 1.795; the M49 densities in Table |[ 
by 2.743; and the NGC 1399 numbers, by 1.875. Once this is done, the total number densities are multiplied by a mean 
globular cluster mass to obtain the mass densities T,ci{Rgc)- 

The mean cluster mass adopted here is (to)c1 = 2.4 x 10^ Mq, which (as was mentioned in §2.1 above) is the value for 
the Milky Way GCS. In both its shape and its mass scale, the Galactic GCLF, or mass function, bears enough similarity 
to those of the cluster systems being discussed here (and, indeed, to most others: Harris 1991) that the same (m)ci may 
safely be applied to them as well. The only issue is whether or not there is any systematic radial trend in mean globular 
cluster masses. It was also pointed out in §2.1 that there are no strong radial dependences observed in the GCLF of M87 
(McLaughlin et al. 1994; Harris et al. 1998), suggesting that (m)ci is constant with rgc in this galaxy, and thus likely 
also in similar systems like M49 and NGC 1399. The question can be addressed in more detail only for the Milky Way, 
which — although clearly an imperfect substitute for these giant ellipticals — is the only galaxy where accurate luminosities 
and three-dimensional galactocentric distances both can be assigned to globulars on an individual basis. 

Figure || shows that there is no dependence on radius in the observed (m)ci of the Milky Way. The top panel plots the 
volume number density of globulars in the Galactic halo, derived from the database of Harris (1996); drawn as the solid 
line is a least-squares fit of a parametric density distribution in the family discussed by Dehnen (1993) and Tremaine et 
al. (1994). (Only metal-poor clusters have been counted, as these are the ones that display true halo kinematics: Zinn 
1985; Minniti 1995.) In the bottom panel, the total V-band luminosities of the clusters (taken again from Harris 1996) 
have been summed separately in each radial bin, and converted to total cluster masses by applying a mass-to-light ratio of 
Tv ~ 2 (Mandushev, Spassova, & Staneva 1991; Pryor & Meylan 1993). The mass densities pc\ that result are shown as 
the circular points. The solid line now corresponds to the same nd fit from the top panel, but multiplied by the average 
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Table 2 

Density Profile of the M87 GCS^ 



Radius 


Rgc 






nci 


(arcmin) 


(avc.) 


(arcmin^^) 


(ave.) 


(arcmin""^) 


0.153 - 


0.328 


0.244 


282.6±55.8 


0.245 


243.7±155.6 


0.328 - 


0.655 


0.494 


197.2±19.9 


0.505 


89.34±34.01 


0.655 


0.983 


0.808 


145.3±16.5 


0.827 


68.58±20.55 


0.983 


1.310 


1.157 


89.10±12.40 


1.152 


18.16±12.94 


1.310 


1.870 


1.500 


74.46±8.78 


1.602 


19.33±5.23 


1.870 - 


2.160 


2.010 


48.92±4.25 


2.018 


11.36±3.56 


2.160 - 


2.490 


2.320 


37.81±3.26 


2.328 


4.015±2.391 


2.490 - 


2.880 


2.680 


35.75±2.80 


2.689 


6.404±1.722 


2.880 - 


3.320 


3.090 


26.93±2.05 


3.104 


3.423±1.183 


3.320 - 


3.830 


3.570 


21.74±1.65 


3.580 


2.875±0.818 


3.830 - 


4.430 


4.120 


16.78±1.34 


4.135 


0.926±0.587 


4.430 - 


5.110 


4.760 


16.38±1.23 


4.776 


1.447±0.469 


5.110 - 


5.900 


5.490 


13.26±1.13 


5.512 


1.202±0.360 


5.900 - 


6.820 


6.350 


9.79±0.93 


6.368 


0.768±0.271 


6.820 - 


7.870 


7.330 


7.16±0.96 


7.354 


0.442±0.219 


7.870 - 


8.973 


8 403 


K ao+n 

iJ.OZi_l_\J. i fj 


8 431 




8.973 - 


9.962 


9.454 


5.68±0.69 


9.474 


0.346±0.154 


9.962 - 


10.95 


10.44 


4.20±0.65 


10.46 


0.222±0.135 


10.95 - 


11.93 


11.43 


3.32±0.61 


11.45 


0.079±0.124 


11.93 - 


12.92 


12.42 


3.23±0.61 


12.43 


0.156±0.120 


12.92 - 


13.90 


13.40 


2.55±0.58 


13.41 


0.008±0.113 


13.90 - 


14.87 


14.38 


3.05±0.57 


14.39 


0.067±0.104 


14.87- 


15.84 


15.35 


3.27±0.57 


15.36 


0.143±0.103 


15.84 - 


16.81 


16.32 


2.59±0.55 


16.33 


0.212±0.093 


16.81 - 


22.52 


19.46 


1.01±0.43 


19.77 


0.024±0.010 



^Limiting magnitude Viim = 23.9; divide by 0.5572 to obtain 
the density over all magnitudes. 



Table 3 

Density Profile of the M49 GCS^ 



Radius 


i?gc 














(arcmin) 


(avc.) 


(arcmin ^) 


(avc.) 


( 


arcmin 




0.732 - 


1.464 


1.035 


29 


.83±2.53 


1.129 


8, 


,041±1, 


.684 


1.464 - 


2.196 


1.793 


18 


.00±1.61 


1.848 


3 


,818±0, 


.874 


2.196 - 


2.928 


2.536 


11 


.40±1.11 


2.575 


1, 


.163±0, 


.543 


2.928 - 


3.660 


3.274 


9. 


,82±1.06 


3.304 


1, 


.307±0, 


.426 


3.660 - 


4.392 


4.009 


7. 


,02±0.80 


4.034 


0, 


.759±0, 


,292 


4.392 ~ 


4.996 


4.684 


5. 


,48±0.63 


4.699 


0, 


.464±0, 


.279 


4.996 - 


5.992 


5.471 


4. 


,51±0.59 


5.505 


0, 


.384±0, 


.169 


5.992 - 


6.988 


6.471 


3. 


,29±0.52 


6.500 


0, 


.276±0, 


.138 


6.988 - 


7.981 


7.467 


2. 


,30±0.48 


7.493 


0, 


.120±0, 


.120 


7.981 - 


8.973 


8.462 


1. 


,97±0.45 


8.484 


0, 


.092±0, 


.097 


8.973 - 


10.95 


9.912 


1. 


,62±0.31 


9.985 


0, 


.070±0, 


.044 


10.95 ~ 


13.90 


12.33 


1. 


,18±0.23 


12.47 


0, 


.054±0, 


.026 


13.90 - 


16.81 


15.28 


0. 


,66±0.22 


15.39 


0, 


.028±0, 


.023 


16.81 - 


19.69 


18.19 


0. 


,37±0.21 


18.28 


0, 


.010±0, 


.005 



^Limiting magnitude Mim = 23.3; divide by 0.3646 to obtain 
the density over all magnitudes. 
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mass calculated for all globulars with 2kpc < Tgc < 40kpc. This shows that pci = ("i)ci"-ci to a good approximation at 
every r^^i a-nd therefore that the mean globular cluster mass is not a strong function of galactocentric radius. The other 
curves in the bottom panel of this Figure are discussed in §4.3. For now, the point is simply to lend some circumstantial 
support to the claim that (m)ci ~ 2.4 x 10^ Mq should not vary significantly with galactocentric radius in M87, M49, and 
NGC 1399. 

The total projected mass densities for these three GCSs are shown as the filled circles in Fig. ||, where angular distances 
and areas have been converted to linear units using Z3(Virgo) — 15 Mpc and Z)(Fornax) = 16.5 Mpc. The final Ed have 
also been scaled upwards by additional factors of 300 — 400 for this Figure, to facilitate direct comparisons with the mass 
surface densities of the galaxies' stars (drawn as lines). These latter distributions are derived from published i?-band 
surface brightness profiles (M87: de Vaucouleurs & Nieto 1978, 1979; M49: King 1978, Caon et al. 1994; NGC 1399: Caon 
et al. 1994), according to the standard relation 

Sstars/(MQ pC-2) = Ts X IQ^M27 .OB+Au-^.u) (12) 

for /is in units of mag arcsec"^. (Extinctions of Ab = 0.09 mag, 0, and are adopted for M87, M49, and NGC 1399; 
see Burstein & Heiles 1984.) The stellar mass-to-light ratios for all three galaxies come from van der Marel's (1991) 
kinematic analysis of their cores; his values have been corrected for the distances D assumed here, and are taken to be 
independent of Rgc- 

The comparisons of Sci and Sgtars in Fig. ^ illustrate the Sn "problems" in M87 and show that, contrary to previous 
claims, they do not exist in M49 and NGC 1399. 

As was indicated above, the global specific frequency of M49 is quite normal; the most recent estimate places it at 
Sn = 4.7 ± 0.6 (Lee et al. 1998). Given a stellar mass-to-light ratio of Ty = 7.5 [which follows from Tb — 10.4 and 
{B — V) — 1.0], equation (||) then implies that Mgcs/Mgtars — (1.7 ± 0.2) x 10"'^ globally in this galaxy. In keeping with 
this, the middle panel of Fig. shows directly that Sci/Estars — 2.3 x 10~^ at essentially any Rgc ^ 15 kpc. This is 
slightly higher than the globalTJCS mass fraction because the globular cluster surface density profile is shallower than 
the stellar distribution at small galactocentric radii. As was discussed in §2.2, such a feature might result, at least in 
part, from the dynamical destruction of globulars in the densest regions of the galaxy. [Again, the theoretical work of 
Murali & Weinberg 1997b, or of Aguilar et al. 1988, suggests that the important spatial scale is generally of order a 
stellar effective radius, i?eff . These are marked by vertical lines for the galaxies in Fig. |^; note that i?eff(M49) ~ 16 kpc] 
This aside, M49 has been put forward as an example of the second specific frequency problem (beginning with Harris 
1986, and continuing on to Lee et al. 1998) on the basis of comparisons at large Rgc ^ ^cff between Eci and the Sgtars 
derived from the surface photometry of King (1978). King's data are shown as the dashed line in the middle of Fig. ^ 
and they clearly do suggest that the projected stellar density falls off more steeply with radius than the GCS density in 
M49. However, the solid line in this panel represents the more recent /i^ data of Caon et al. (1994), which include many 
more measurements at Rgc ^ 30 kpc and do not show any GCS-halo discrepancy. The local specific frequency, which 



Table 4 

Density Profile of the NGC 1399 GCS^ 





Radius 


i?gc 




Ncx 




Tgc 










(arcmin) 


(ave.) 


(arcmin 




(ave.) 


( 


arcmin 




0, 


.761 ~ 


1.518 


1.075 


54 


.88±4 


.60 


1.171 


13.74±0 


.31 


1, 


.518 - 


2.274 


1.858 


35 


.30±3 


.10 


1.915 


7, 


.719±0, 


,163 


2. 


.274- 


3.031 


2.625 


20 


.64±2 


.11 


2.666 


4, 


.288±0, 


,949 


3, 


.031 - 


3.788 


3.388 


10 


.81±1 


.48 


3.420 


1, 


.141±0, 


,607 


3, 


.788 - 


4.544 


4.149 


8. 


,77±1. 


.27 


4.175 


0, 


.957±0, 


,472 


4, 


.544- 


5.301 


4.908 


6. 


,72±1. 


.10 


4.930 


0, 


.618±0, 


,393 


5, 


.301 - 


6.058 


5.667 


5. 


,31±0. 


.97 


5.686 


0, 


.496±0, 


,314 


6, 


.058 - 


6.814 


6.425 


3. 


,93±0. 


.86 


6.442 


0, 


.402±0, 


,267 


6, 


.814- 


7.571 


7.183 


2. 


,61±0. 


,79 


7.197 


0, 


.108±0, 


,261 


7, 


.571 - 


8.104 


7.833 


2. 


,60±0. 


,99 


7.840 


0, 


.093±0, 


,332 


8, 


.104- 


9.573 


8.808 


2. 


,32±0. 


,71 


8.854 


0, 


.175±0, 


,134 


9, 


.573 - 


11.04 


10.21 


1. 


,39±0. 


,68 


10.32 


0, 


.045±0, 


,111 


11.04- 


12.50 


11.75 


1. 


,32±0. 


,55 


11.78 


0, 


.095±0, 


,097 


12.50 - 


13.97 


13.21 


0. 


,65±0. 


,62 


13.24 


0, 


.024±0, 


,023 



^Limiting magnitude Vnni — 24.0; divide by 0.5332 to obtain 
the density over all magnitudes. 
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r/kpc 

Fig. 2. — Volume number- and mass-density profiles of the metal-poor (halo) globular clusters in the Milky Way. Individual magnitudes, 
metallicities, and galactocentric radii for 112 clusters are taken from the compilation of Harris (1996). Solid curves represent a parametric fit 
to the n^i distribution, which is multiplied by a mean cluster mass of {m)^i = 2.4 X IQ^Mq to obtain the mass densities p^. The good fit to 
the data in both the top and bottom panels shows that (m)ci is not a function of Tgc in the Milky Way. Broken curves in the bottom panel 
refer to the density profile Pstars(?'gc) of the Pop. II spheroid; see §4.3 of the text. 

is proportional to Eci/5^stars, then takes on a roughly constant value throughout the outer halo of M49, and the second 
Sn problem does not appear there.|^ The differences between the light profiles of King (1978) and Caon et al. (1994) 
illustrate the difficulty of accurate sky subtraction at faint surface brightness levels, and make the point that background 
corrections are an important consideration in any discussion of galaxy (or GCS) structure on large spatial scales. The 
surface photometry of Caon et al. (1994) is adopted here for M49, and thus Ed oc Egtars at large Rgc, but this important 
caveat must be kept in mind. 

In §3.2, it is shown that the gas in M49 adds up to a small fraction of its stellar mass, in which case the local cluster 
formation efficiency can be estimated by tci — Pci/(Pgas + Pstars) — Pci/Pstars- The fact that Sd oc Sgtars therefore implies 
that Ed was essentially independent of galactocentric radius beyond about i?eff — 16 kpc. An extrapolation of this 
result to smaller distances requires that the systematic decrease in Eci/Estars there has resulted entirely from the gradual 
destruction of (initially bound) giobulars over a Hubble time. As possible anecdotal support for this interpretation, recall 
that globally Mgcs/Mgtars — 0.0017 for M49, while locally Sd/Egtars — 0.0023 at large radii. If the latter number is 
assumed to have held all the way to the center of the galaxy initially, then a global Mg"g*/Afg"^Jg ~ 0.0023 is implied. In 
such a scenario, if Mgtars were roughly conserved over a Hubble time while giobulars were dynamically disrupted (their 
remains adding nothing significant to the total stellar mass), the global ratio of current to initial GCS masses would be 
Mgcs/M™g* ~ 0.0017/0.0023 ~ 0.75. This is surprisingly close to the result of the general GCS mass-function arguments 
in §2.1 (see eq. S). However, this could be just a coincidence; it remains possible that ed really was lower at small 
galactocentric radii, i.e., that the GCS formed with a "core" which was somewhat larger than that of the unclustered 
stellar distribution, and which was only enhanced over time. The current data cannot settle the question either way. 

The bottom panel of Fig. ^ shows that the second Sn problem does not exist beyond a stellar i?ofr in NGC 1399, either. 
Every study of this GCS has appreciated this point, and it is not a new result (see the references cited above). (Note 
also that the surface brightness profile used here for NGC 1399, from Caon et al. 1994, is consistent with that of Killeen 
& Bicknell 1988 at large -Rgc-) On the other hand, the global Sn in this galaxy — and therefore the local Ed/Estais at 

^Lee et al. (1998) have recently claimed evidence for a difference in the radial distributions Ed of red (metal-rich) and blue (or metal-poor) 
globular clusters in M49. They argue that the former are more centrally concentrated, i.e., exhibit a steeper fall-off in density with Rgc, than 
the latter. They also suggest that the "red" component of the M49 GCS follows the galaxy's stellar profile Sstars, while the "blue" component 
is spatially more extended (see also Forbes et al. 1997). However, since Lee et al. were unable to determine the background densities Nf, 
separately in the blue and red subsets that they defined for the GCS, further observations are required to verify their claim (and comparisons 
then should perhaps be made with the surface brightness profile of Caon et al. 1994). Moreover, any real differences in the volume density 
profiles of red and blue giobulars will be artificially amplified in projection. If Lee et al.'s result is confirmed, it will be important to understand 
in detail; but allot the GCSs discussed here are considered only in their entireties, without any attempt to distinguish between subpopulations. 
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Fig. 3. — Comparison of the projected mass density profiles of globular clusters (points) with those of the unclustered halo stars (lines) 
in M87, M49, and NGC 1399. GCS densities have been corrected for photometric incompleteness (see text), and scaled up by the inverse of 
the cluster formation efficiencies derived in §3.2 below. Mass-to-light ratios for the stars are taken in each galaxy from van der Marel (1991). 
The broken vertical lines mark the effective radius of each galaxy: i?eff = 7.0 kpc for M87 (de Vaucouleurs & Nieto 1978), 16.4 kpc for M49 
(Caon et al. 1994), and 10.2 kpc in NGC 1399 (Caon et al. 1994). The disparity between Sd and Sstars in M87 (top panel) embodies the 
specific frequency "problems" there. 

any single i?gc — has always been thought to be significantly higher than the normal, M49-like value of ^5: the latest 
calculation claims Sn = 12 ± 3 (Kissler-Patig et al. 1997).^ However, it is immediately apparent from Fig. ^ that this is 
not the case: the bulk scaling of 1/0.0029 required to match the GCS radial distribution to the stellar mass profile in NGC 
1399 is identical, within the observational uncertainties, to the 1/0.0023 needed in M49. Since the gas mass is negligible 
compared to the stellar mass on < 70-kpc scales in NGC 1399, as it is in M49, EcI — Sd/Sstars holds for both systems, and 
their cluster formation efficiencies are inferred to have been essentially identical. The spuriously high global Sn values 
typically quoted for NGC 1399 can be traced back to the practice of assigning an absolute magnitude of My — —21.6 to 
the galaxy (from the Third Reference Catalogue: de Vaucouleurs et al. 1991); but, in fact. My ~ —22.1 is indicated by 
direct integration (to Rgc = 70 kpc, assuming D — 16.5 Mpc) of the surface brightness profiles of either Caon et al. (1994) 
or Killeen & Bicknell (1988). Given this brighter magnitude, and an estimated GCS population of A/tot = 4 700± 400 
projected to within 70 kpc of NGC 1399 (see §5), equation (|l|) yields S'Ar(NGC 1399) = 7.0 ± 0.6, much closer to the 
value for M49 (see also Ostrov et al. 1998). This is consistent with the alternate representation of the data in Fig. |^, and 
confirms that the first specific frequency problem is not an issue in NGC 1399. 

It appears, therefore, that there are no Sn problems in either M49 or NGC 1399. There is no evidence for any 
fundamental variation in ed from the total populations of these two GCSs, and the local Cci in each is independent of 
galactocentric radius beyond ~ 1 Rcs- These results lend some weight to the claim (§2.2) that the first and second specific 
frequency problems seem likely to appear, and should be solved, only together. 

Both of these problems are genuine in M87, as the top panel of Fig. ^ attests. (The surface photometry of de Vaucouleurs 
& Nieto 1978, which is the standard, is well sampled and is supported by the independent studies of Oemler 1976 and 
King 1978; and the total My used to calculate the galaxy's Sn is taken from integration of its fiB profile.) The second 
specific frequency problem is evidenced by the shallower decline of Sci compared to Sgtars at Rgc ^ ReS — 7 kpc. This 
also causes the first problem, since scaling the total GCS density upwards by the factors appropriate for M49 or NGC 
1399 obviously puts Sd well in excess of Sgtars throughout most of the galaxy. Indeed, the most recent estimate for the 
global Sn of M87 is 14.1 ± 1.6 (Harris et al. 1998), which is fully a factor of three higher than that of M49. Here it is 
necessary to consider the distribution of hot gas — which is known to exist in great quantities around M87 — in evaluating 
ed from equation (|l0|). 

^After this paper was submitted, a new study of NGC 1399 by Ostrov, Geisler, &: Forte (1998) found that the specific frequency of its GCS 
is a much more modest ^jv = 5.6 it 1.0. This result is consistent with that obtained here, and in fact is even closer to the value for M49; it 
supports the conclusion that NGC 1399 does not exhibit the first Sjv problem. 
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Fig. 4. — Comparison of the (mass) volume density profiles of the GCS (filled circles), halo stars (thin solid line), and X-ray gas (open 
squares and broken line) in M87. The volume densities of globular clusters and field stars follow from deprojections of the surfa ce d er isitv 
profiles in Fig. ^. To deproject the galaxy surface photometry, it is assumed that pstars oc r'^^ at large radii rgc > 100 kpc (see eqs. [ |A4| , Ja^ , 
and ]A8| of the Appendix). The vertical line segment is placed at the projected effective radius of the stellar light, as in Fig. ^. Thick black 
line is tne sum of the stellar and gas densities. Note that pcl is directly proportional to (pstars + Pgas) beyond rgc ~ i?eff> giving a constant 
estimated cluster formation efficiency of t^i = 0.0032 it 0.0005 there. The sharp feature at rgc ~ 20 kpc marks the onset of the cD envelope 
in the stars and GCS of M87; it is visible as a smoother "hump" in projection (de Vaucouleurs & Nieto 1978; McLaughlin et al. 1993). 



3.2. Volume Densities 

The volume density profiles of circumgalactic and intracluster gas are routinely derived in the course of X-ray studies 
of elliptical galaxies. To include such a halo component in this discussion, it is first necessary either to deproject observed 
GCS and stellar surface densities, or to compute the gas surface density via a standard projection integral. In general, 
the first of these options potentially has the most insight to offer, since projected quantities always mix information from 
a wide range of physical radii. 

An Appendix therefore develops a nonparametric deprojection algorithm which is explicitly designed to handle GCS 
data that have been binned in projected radius, Rgc (i-e., counted in circular annuli). This procedure returns a coarse 
volume density profile — Tid vs. three-dimensional rgc — expressed as a series of average densities in a number of concentric 
spherical shells. The algorithm is preferably applied to surface density profiles that extend to fairly large radii (so that 
corrections for cluster populations beyond the observed field of view can be estimated with confidence, and are as small 
as possible), and it makes most sense for A^ci distributions that are monotonically decreasing functions of i?gc (otherwise, 
the derived rid can be zero or negative in some shells) . After re-binning to smooth some of the original data, as described 



above, the composite GCS radial distributions in Tables g through ^ satisfy each of these criteria. Equation (A5) of 
the Appendix has therefore been applied to each of these Nc\ profiles, with f{Ri^i,Ri) (a sort of "residual background" 
correction) calculated by assuming that rid cx beyond the last observed radius of each GCS. (This reflects the fact 
that the projected A'd shows roughly an R~^ behavior at large Rgc in each case. In any event, only the outermost one or 
two points in each deprojectedprofile are affected — and then not greatly — by this choice.) The resulting nd are listed in 
the fifth columns of Tables ||, ^ and ||. The radii in the first columns, which delimit projected annuli in reference to the 
surface density data, are identified with the inner and outer three-dimensional radii of spherical shells in reference to the 
volume densities. 

The deprojection algorithm is an iterative one: the volume density in a spherical shell with some mean radius is 
computed by first subtracting, from the projected density at that radius, the contribution of any globulars that actually 
lie at larger physical rgc. Thus, the uncertainties of the derived rid at any two radii in one system are not independent. 
The (1-cr) errorbars quoted in Tables |[ ||, and |^ are the standard deviations of results from 1000 different deprojections 
in each case. Each trial deprojection in each GCS begins with a different A'd profile, in which the surface density for any 
given annulus is drawn at random from a Gaussian distribution that is centered on the value given in Column 3 of the 
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r/kpc 



Fig. 5. — Same as Fig. W, but for M49. Two separate stellar density profiles are shown in this case, derived by deprojecting the surface 
brightness profiles of King (T_978; dashed line) and of Caon et al. (1994; thin solid line). The projected half-light radius of the latter is marked 
with a vertical line. For both deprojections, it is assumed that pstars oc r'^^ 3't large radii. The bold line is the sum of pstars (from the profile 
of Caon et al.) and Pgaa (open squares and dotted line). At radii r > R^g, the ratio Pcl/(Pstars + Pgas) is again constant with radius, giving 

a best estimate of e^i = 0.0023 it 0.0005. This is within 30% of the value found in M87, and the two are identical within their errors. 



appropriate Table, and that has a dispersion equal to the observational uncertainty there. 

The tabulated number densities Uci have the same limiting magnitudes as the corresponding surface densities, and are 
therefore multiplied by the same two factors (so effectively scaling to Viim = oo, and applying the mean globular cluster 
mass (m)ci — 2.4 x 10^ A/0) to find the total GCS mass densities PcKj^gc)- These distributions are plotted as the filled 
circles in Figs. ^, ^, and ^, where the density in ea ch spherical shell has been put at a radius Tgc, given in Columns 4 of 
Tables ||-^, calculated according to equation ( A10| ). As in Fig. ||, the pci profiles have been scaled up by factors of a few 



hundred in order to compare them directly with the stellar mass densities Pstars (''gc ) (thin solid lines in Figs. ^ to ^). These 
have been obtained by applying the alg orithm of the Appendix to the projected profiles Sstars(^gc) derived above from 
i?-band surface photometry (see eq. |l2|]).p| It is immediately apparent that the good agreement seen in Fig. ^between the 
GCS and stellar densities in NGC 1399 and M49 (for the data of Caon et al. 1994) persists after deprojection, as should 
be the case. More noteworthy is the fact that the discrepancy between pci and Pstars in M87, while certainly still present, 
is somewhat less pronounced than that between Sci and Sstars- For example, the truly local picture given by Fig. H shows 
that the second Sm problem here is confined to galactocentric radii rgc ^ 30-40 kpc, and does not extend inward to the 
~ 15 kpc scales suggested by Fig. ||. This just reflects the fact that the projection of any two volume densities inevitably 
amplifies any small differences that may actually exist between them, and artificially associates disparities at large rgc 
with smaller projected radii. A similar comment applies to the comparison of pci with pstars in M49: it is now clear, 
from Fig. |^, that any argument for the second Sm problem in this galaxy could only be made for three-dimensional radii 
greater than about 50 kpc, and then (as in §3.1) only if King's (1978) ^,3 profile — which includes just a single data point 
at such large rgc — were preferred over that of Caon et al. (1994). 

The gas densities in these systems have been obtained from published electron number densities by way of the relation 

= 2.865 x 10-^ (J^)(-^^) , (13) 



M0PC-3 ■ VO.6O/ VlO-3cm-3 

where a mean particle mass ji ~ 0.6 (in units of uih) is equally appropriate for hot {T ^ 10^-10® K) and highly ionized 

^AU of the surface photometry used here has been averaged in circular annuli, as is required by the deprojection method in the Appendix. 
[Effective circular radii Rgc = (afe)^/^ are assigned to elliptical isophotes with semi-major and -minor axes a and b.] Note that these data are 
generally published just as tables of /ig vs. Hgc, with no information given on the inner and outer radii of any annuli used in the measurements. 



Thus, to apply equation (A5), annuli have been defined in an ad hoc manner, by setting boundaries at the geometric means (filfi2)^/2 of 
successive pairs of projected radii in the published data. Data points are generally spaced so closely as to make the difference between this and 
any other choice unimportant. 
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Fig. 6. — Same as Figs. yand|5[ but for NGC 1399. Filled circles are again the GCS densities; thin solid line, the stellar densities (obtained 
by deprojection of B-band surface brightnesses, correcting for pstars oc r'^^ beyond about 70 kpc); dotted line, the gas density; and bold solid 
line, the sum (pgas + Pstars)- As before, the vertical line is placed at the stellar i?cff- The estimated cluster formation efficiency is again 
independent of galactocentric radius: t^i = 0.0029 ± 0.0008. 



plasmas of any composition from primordial to solar. 

The open squares for the M87 gas in Fig. ^ come from scaling the electron densities of Nulsen & Bohringer (1995) for 
the different Virgo distance assumed here (15 Mpc vs. their 20), according to oc £1-1/2. (These data were published as 
volume densities after Nulsen & Bohringer applied a deprojection algorithm, similar to that derived in the Appendix, to 
their i?05^r spectra.) The dotted line in Fig. | is a fit taken from McLaughhn (1999): 



Pgas(M87) = 5.60 X 10"^ Mq pc" 



1070 kpc 



1070 kpc 



(14) 



so that, essentially, pgas c>c r"^.^ over the area of interest here. This fit has been used to compute the total pgas + Pstars at 
every radius with a measured stellar density; that sum is drawn as the bold solid line in the Figure. 

Similarly, the open squares in Fig. here use the Tig read from Irwin & Sarazin's (1996) Figure 6, after scaling to 
D — \b Mpc from their assumed M49 distance of 25.8 Mpc. The dotted line in this Figure is a scaled fit (from Brighenti 
& Mathews 1998) to these data plus those of Trinchieri et al. (1986): 



Pgas(M49) 



3.57 X 10 



-3 



2.24 X 10" 



1.50 X 10" 



1 + (rgc/0.807kpc)2 1 + (rgc/7.18 kpc) 



1.14 



l + (rgc/75.6kpc) 



1.19 



Mqpc- 



(15) 



The heavy black line again represents the sum of pgas + Pstars- 

The broken line for pgas in NGC 1399 (Fig. ph comes from Tsai's (1993) fit to the X-ray observations of Killeen & 
Bicknell (1988), corrected for I?(Fornax) = 16.5 Mpc: 



Pgas(NGC 1399) = 1.03 x IO^^Mqpc 



-3 



' gc 



0.304 kpc 



-0.615 



(16) 



It is clear that stars outweigh the gas everywhere in M49 and NGC 1399: locally, Pgas/pstars <, 0.3 even at the largest 
radii in Figs. ^ and ||; globally, Mgas/Mgtars ^ 0.05 within rgc = 100 kpc and 70 kpc respectively. In both systems, then. 



Mgcs/Mstars and (albeit to shghtly lower accuracy at very large rgc) Pci/ (pgas + Pstars) — Pci/Ps 



It is indeed valid to use standard specific frequencies as direct measures of globular cluster formation efficiencies in these 
two cases, as was done at the end of §3.1 above. To be complete, however, all of the stars, gas, and globular clusters can 
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now be used to estimate ed more precisely, according to the basic equation (|1C|). By summing the data at large Tgc in 
each of M49 and NGC 1399, then, average Cci values are obtained: 

Mrrcsirzc > 12.5 kpc) , , , , 
s^^^ ^ = 0.0023 ± 0.0005 (M49) (17) 



Mga,s{rgc > 12.5kpc) + Mstars(?'gc > 12.5 kpc) 
and 

Mgcs(?^gc > 11 kpc) 



,x . N ,r / -= 0.0029 ±0.0008 (NGC 1399), (18) 

Mgas(rgc > llkpc) + A/,tars(?-gc > Ukpc) 

which account for the scalings applied to Eci and pci in Figs. ^, ^, and ^. As has already been discussed, these are equally 
well interpreted as either local or global formation efficiencies, since the ratio Pci/(Pgas + Pstars) shows no variations with 
Tgc beyond a stellar effective radius (or even inside this, in NGC 1399), i.e., outside those regions of the GCSs which could 
be significantly dynamically evolved. (As in Fig. 0, the stellar R^g are marked with vertical broken lines in Figs. B, 0, and 

I) n 

The situation is somewhat different for M87: referring to Fig. |4|, Pgas/pstars > 1 at all radii Tgc > 40 kpc (which is where 

the second Sn problem begins in the deprojected data), and Mgas/Mgtars — 0.35-0.40 globally (i.e., integrated over scales 
Tgc < 100 kpc). The hot gas around M87 is therefore an important factor in calculating its cluster formation efficiency. 
As before, outside of the central regions of the galaxy it is found that 



A'/gcs(fgc > 11 kpc) 



A^gas(?'gc > 11 kpc) + Mstars('^gc > H kpc) 



0.0032 ±0.0005 (M87) , (19) 



explaining the scaling factor applied to the M87 GCS densities in Figs. ^ and ^. Figure H shows that this average 
value is valid also as the local Cci = Pci/{pgas ± Pstars) at any galactocentric position beyond about ReS = 7 kpc; as 
in M49 and NGC 1399, the inferred globular cluster formation efficiency in M87 was independent of rgc at large radii. 
In this case, however, if the gas were not included in the calculation, much larger — and radially varying — formation 
efficiencies would have been indicated: for example, Mgcs(?'gc > 11 kpc)/Mstars(?'gc > 11 kpc) — 0.0056 ± 0.0008, and 
Mgcs(rgc > 40kpc)/Afstars(?'gc > 40 kpc) = 0.0095 ± 0.0023. Note that the first of these ratios, A/gcs/M^tars = 0.0056 
for Tgc > 11 kpc, is a factor of 2.4 higher than the corresponding quantity in M49. In combination with the different 
stellar mass-to-light ratios in these two galaxies (Tg = 14.6 vs. 10.4; see van der Marel 1991, or Fig. || above), this just 
corresponds to the factor of 3 difference observed in their global Sn values. The high local and global specific frequencies 
in M87 are thus seen to result from the presence of a population of globular clusters that are intimately associated with 
the hot gas there. 

3.3. Gas Surface Densities 

As was mentioned earlier, an alternate comparison of a galaxy's gas distribution with its stars and globular clusters could 
make use of only surface densities — Estars, ^ci, and a derived Sgas — as functions of projected Rgc- The advantage of such 
an approach is that Egas is on an equal physical footing with Sgtars and Ed, and the latter are basically unmanipulated 
data with smaller relative uncertainties than the derived Pstars and pd profiles. The cost, aside from possible uncertainties 
in any extrapolations required to obtain Egas from Pgas, is that differences between projected distributions are not easily 
related to local effects at some true, three-dimensional radius (e.g., recall the comparison between projected Ed vs. Egtars 
and deprojected pd vs. Pstais in M87). More to the point, if the ratio Ed/(Egas + Egtars) varies as a function of Rgc in 
a galaxy, then variations in the physical quantity of real interest, the local ed = Pd/(Pgas + Pstars), are clearly indicated, 
but cannot be properly quantified without resorting to a deprojection of some sort. For this reason, the procedure of §3.2 
must be favored as the best way, generally speaking, to measure local ed. 

That said, however, in the galaxies studied here, ed is found to be independent of radius at large rgc, i.e., pd oc 
(Pgas + Pstais)- In this specific instance, it should then also happen that Ed cx (Egas + Egtars) at large i?gc, with a constant 
of proportionality that is exactly Cd. It is therefore useful to compute Egas and the projected analogue of Cci in M87, 
M49, and NGC 1399, in order both to verify the results of §3.2 (which depend in good part on the behavior of pc\ at large 
Tgc, where the deprojected data are noisiest) and to obtain potentially higher-precision estimates of the cluster formation 
efficiencies in these systems. 



First, then, the projected mass densities of gas are obtained from the volume density profiles in equations (14), ( [15D 
and (E^, by way of the standard integral 



Egas(-Rgc) = 2 / Pgas{rgc)dz , (20) 



where the coordinate z — {r"^^ — Rg^Y^^ measures distance along the line of sight, and r^ax represents the total physical 
extent of the X-ray corona. For M87, this is taken to be the virial radius of the Virgo Cluster: Tmax — 1-5 Mpc (see 
McLaughlin 1999). For M49, which is not at the center of Virgo, r^nax = 500 kpc is adopted. NGC 1399 lies at the center 
of the Fornax Cluster, which is somewhat less massive than Virgo (e.g., Jones et al. 1997), and rmax is set at 1 Mpc there. 
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R/kpc 

Fig. 7. — Comparison of GCS, stellar, and gas surface densities in MS'L. M49, and NGC ]|399- Xiie bold solid line in each panel traces the 
sum Egas + Satars, and the galaxies' effective radii are marked as in Fig. pi Also as in Figs. Hand Q the solid and dashed lines for the stellar 
densities in M49 refer to the surface photometry of Caon et al. (1994) and King (1978), respectively. The implied values of are essentially 
the same as those derived from the volume density profiles in all three galaxies, although the ~ 20% lower value in M87 may be significant. 
The cluster formation efficiency is again inferred to have been locally invariant in each of these systems (at least beyond i?gc ~ iJcff ) i a^nd not 
significantly different from one to the other. 



There is a certain arbitrariness to each of these choices for Tmax, but any of them can be changed by a factor of a few 
without affecting the conclusions which follow. 

The resulting Egas profiles are shown in Fig. along with the same Sstars and Sci profiles plotted in Fig. ^ above 
(although the Sci have been scaled up by slightly different factors in the current Figure). It is easily seen that Sci esc 
(Egas + Egtars) at large Rgc in each galaxy, just as expected following the analysis in §3.2. Moreover, the constants of 
proportionality are in reasonable agreement with the Cci derived from the deprojected data: by taking the direct mean of 
the surface-density ratios beyond the galaxies' effective radii (so as once again to avoid a zone of potential GCS erosion), 
it is found that 

V. V r 0.0026 ±0.0003 , M87 

) = <^ 0.0022 ± 0.0004 , M49 (21) 

Sgas + Sstars / _R^,^fl,„ [ 0.0026 ± 0.0006 , NGC 1399 . 

The result for M49 is indistinguishable from that obtained in either §3.1 or §3.2, owing to the truly low gas densities in 
this galaxy. Similarly, the projected Cci in NGC 1399 is identical, within the uncertainties, to the deprojected quantity 
Pci/(Pgas + Pstars) = 0.0029 ± 0.0008. The number in equation ( |2l] ) is formally 10% smaller because of the apparently large 
projected gas densities at radii Rgc ^ 40-50 kpc in NGC 1399. However, these values are almost certainly overestimates: 
in deriving Egas, the shallow volume density profile of equation (|l^), which roughly has Pgas cx: r^c''^, has been assumed 
to hold even on Mpc spatial scales in Fornax — well beyond the range of actual X-ray observations — and this is likely 
unreaHstic. On the other hand, the projected ?ci = 0.0026 ±0.0003 indicated for M87 in Fig. ^ might genuinely be sHghtly 
smaller than the 0.0032 ± 0.0005 inferred from deprojected quantities. The 20% difference between these two numbers is 
significant at just the l-cr level, however, and — although it is interesting that the former is actually closer to the values 
for M49 and NGC 1399 — it is not obvious which is more reliable. The X-ray gas density around M87 (eq. ||l^) has 
again been extrapolated to unobserved radii to compute Egas, but in this case the basic profile actually follows the dark 
matter distribution in Virgo (McLaughhn 1999), and exhibits a realistic, steepening slope at large rgc. Thus, although it 
is possible that Egas is overestimated at large Rgc in M87, any argument for this is not as compelling as in NGC 1399. 

The best measure of eci in M87 will therefore be taken as the error- weighted mean of the values (eqs. jl^ and [^ ) 
obtained by analyzing the deprojected and projected data there; in M49, the choice is of no consequence, but for definite- 
ness the deprojected Cci will be adopted; and the deprojected estimate will also be adopted for NGC 1399, to avoid any 
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concerns that its gas surface density might have been overestimated at large -Rgc- Thus, 

r 0.0028 ±0.0004 , M87 
?ci = < 0.0023 ± 0.0005 , M49 (22) 
i 0.0029 ± 0.0008 , NGC 1399 . 

Once again, these can be seen as estimates of either local or global cluster formation efficiencies if the inner several kpc 
of M87 and M49 are excluded from the discussion. 

4. A UNIVERSAL EFFICIENCY FOR CLUSTER FORMATION 

The results of the previous Section clearly show that working in terms of stellar masses rather than luminosities and 
including the gas in M87 removes both the first and second specific frequency problems there: in all three of M87, M49, and 
NGC 1399, it is found that pd cx (pgas + Pstars), in identical proportions within the uncertainties, outside of the central, 
possibly dynamically depleted regions of the GCSs. (At the most, differences between the constants of proportionality in 
M87 and M49 might remain at the ^--^20% level; but again, this is easily taken up in the observational uncertainties, and 
in any case is a significant improvement over the 300% difference in their specific frequencies.) The interpretation put to 
this is that globular clusters formed with an essentially constant efficiency throughout all of these galaxies, at the level of 

M™^' M 

eci(globular) = —1^ ~ = 0.0026 ± 0.0005 , (23) 

M'gTs Mga. + Mstars ^ ' 



which is the error- weighted mean of equation (^2|) . 

It should be kept in mind that these are all large galaxies with complex evolutionary histories. Thus, as was discussed 
in §2.1, what has been measured here is necessarily a mass-weighted average of ed, possibly over multiple distinct episodes 
of in situ cluster formation and/or of the accretion of globulars formed elsewhere. Nevertheless, the uniformity of the 
observational estimates Cci in these systems — both as a function of rgc within each and from one entire system to another — 
suggests that these concerns may not be too severe, and that equation ( |23| ) may still serve as a reliable guide to a single, 
possibly universal efficiency for cluster formation. This impression can now be checked quantitatively, albeit to rather 
lower precision than has proved possible with the excellent datasets for M87, M49, and NGC 1399, by appealing to 
observations of many other GCSs. All indications are that Cci is indeed similar, at least to first order, from galaxy to 
galaxy and from the protogalactic epoch to the present. 

4.1. Other Ellipticals and BCGs 

Given the definition for EcI (eq. [0), and denoting by G the global ratio Afgas/Aistars in a galaxy, the total population 
of any GCS may be written as 



AA,o, = 2.92xlO«?,(l + G) L.,':^.. . (24) 



IQ^^Lq) {tMqLq^} V2.4xlO5Af0 

This can, of course, also be derived from equation (|^) above, which explicitly relates Sn to e~ci when G = 0. 

Now, in the general discussion of §2, the stellar mass-to-light ratio was taken to be roughly similar from galaxy to 
galaxy. This is not exactly true, however: as was seen in §3, for example, the recovery of similar ed in M87 and M49 from 
their significantly different global and local Sn depends on the ~40% larger Ts of M87 as well as on its greater gas mass. 
Thus, if other ellipticals are to be brought properly into this discussion, it is necessary first to account for the fact that T 
depends systematically on galaxy luminosity. This is a well known result of many fundamental-plane analyses of bright 
ellipticals, and it can be incorporated quantitatively in this analysis by referring to the study of van der Marel (1991; and 
references therein). Specifically, for a sample of 37 early- type systems, van der Marel obtains a mean i?-band mass-to- light 
ratio of T^,: — 3.32/i5o Mq Lq^, corresponding to Tb — 5.93/i5o, at a mean magnitude of Mb — —22.22 -I- 5 log ^,50. For a 
galaxy color of {B — V) = 1.0 (typical of large ellipticals) and Hq = 70 km s~^ Mpc~^, this result corresponds to Ty — 6.0 
at an absolute magnitude of My — —22.5, or Ly.gai — 8.5 x 10^° Lq. In addition to this, van der Marel (1991) finds 
that T increases with Lgai as a power law with exponent 0.35 ± 0.05. Since other studies of the fundamental plane have 
claimed a slightly shallower exponent of 0.25-0.30 (e.g., Faber et al. 1987; Pahre, Djorgovski, & de Carvalho 1995, and 
references therein), Ty cx iygai is adopted here. Thus, 

_:^=Q,3(-^X^Y\ (25) 

It should be noted again that these mass-to-light ratios refer to the stell ar p opulations in the cores of the galaxies, and are 
not unduly influenced by any dark matter on larger scales. If equation ( |2^ ) is used in (p4|), the dependence A/tot « iyg^i, 
or Sn (X iygai: emerges naturally for constant Cd and G. This is essentially the scaling found empirically by Kissler-Patig 
(1997). 
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The other point of concern is, of course, the hot gas content of elhpticals in generaL To get a rough handle on this, 
notice that the (soft) X-ray luminosities of early- type systems scale with their blue luminosities as Lx oc Lg, with x ~ 2-3 
(and with significant scatter about the mean trend: e.g., Forman et al. 1985; Brown & Bregman 1998). This scaling can 

be related to the ratio G — Mgas/Mgtars by noting that Lx esc J p'^^^Tgls dV oc M^^^Tgls for the bremsstrahlung 
emission that dominates Lx in this context. Also, the slope x of the Lx-Lb relation is the same if Ly is used instead, so 
the imphcation is that cx T'^J'^T''^ L^y'"^ . To go further, Tgas bears some relation to the stellar velocity dispersion: 
Tgas c!c a^y , with y ~ 0.7-1.5 (Davis & White 1996; Brown & Bregman 1998); the F-band fundamental plane gives for 
the stellar quantities, tr cx L^^^^R'^-'^^ and R cx L\^^ (e.g., de Carvalho & Djorgovski 1989, 1992; Pahre et al. 1995); and, 
from equation (|25| ) just above, Ty cx Ly^. Putting all of this together, G^ cx iJ^^-'^~^-'^^y ^ Evidently, the uncertainties 
in X and y are not crippling to this analysis, as G cx Ly^*°'^^ is fairly representative of the range of possibilities. The 
scaling can be normalized to apply specifically to brightest cluster galaxies (BCGs) by using the observations of M87, 
which lies at the spatial and dynamical center of the Virgo Cluster. q| As was mentioned in §3.2, the global (rgc < 100 
kpc) ratio of gas to stars in M87 is ~ 0.4 by mass. Since the l^-band luminosity of the galaxy is about 8 x IQ^^Lq (from 
de Vaucouleurs & Nieto 1978, for Z) = 15 Mpc), it is then reasonable to set 

Putting this and equation ( p5| ) both into equation (H), and holding eci fixed, it can be seen that the mean trend expected 
for A/tot between Lv,g&\ = 3 x 10^"Lq and 2 x 10^^L^(so between an imphed G of ~0.1 and 1.6) is roughly Ntot oc i^gai- 
Despite the crudeness of the estimate of G used here, the power 1.75 is just that inferred from observational scalings 
presented by Harris et al. (1998) for a large sample of BCGs. 

The full relation between Atot and Ly,gai is drawn as the heavy black line in Fig. ||, where the total CCS populations 
of 65 early-type galaxies are plotted against their V luminosities. These data are the same as those plotted in Figure 11 
of Harris et al. (1998), and the original references are given in that paper. The open circles in Fig. here denote BCGs 
in a large number of rich clusters and a few poor groups; the filled symbols refer to other ellipticals, some in clusters and 
some more isolated. The bold curve has eci = Afgcs/(Afgas + -Mstars) = 0.0026, constant in all galaxies. Also shown, as 



the dashed line, is the result of keeping this uniform eci, but setting G = in equation (24) while still having Ty increase 
with luminosity as in equation (|25|). The three dotted straight lines are loci of constant specific frequency. A/tot oc Ly,gai, 
as given by equation (l|) for (from top to bottom in the Figure) Sjq = 15, 5, and 1.5. 

It is immediately obvious why a roughly constant Sx — 5 is often assumed for large ellipticals such as these. It is 
equally clear, however, that such a simple characterization tends to overpredict A/tot at low Lygai, and systematically 
underestimates it at high luminosities. Accounting for Ty cx ^y gai goes some of the way towards explaining this while 
allowing for a universal formation efficiency Cci; and referring the GCS populations in BCGs to the galaxies' total stellar 
plus gas masses seems indeed to remove the rest of the first specific frequency problem in general. 

It is encouraging that these rather broad considerations lead to a significantly improved description of the observed 
A/tot ~iy,gai dependence, in which there are no obvious systematic deviations from a uniform Cci — 0.0026 (see also §4.4). 
There remains a good deal of scatter about the mean line in Fig. 0, but this could well reflect departures of individual 
galaxies from the fundamental plane, the mean Lx-Lb (or G-Ly gai) relation, or the assumed {B — V) = 1.0, rather than 
any large variations in the fundamental efficiency of cluster formation. Note, for example, that the observed mass-to-light 
ratio of M87 is Ty ~ 10.5, some 80% higher than that expected from equation (|25|). Correcting for this would bring the 
?ci — 0.0026 curve in Fig. ^ just up to the M87 data point there. To rather a lesser degree, some of the scatter could also 
arise from GCS-to-GGS differences in the relative magnitude of any dynamical reduction of Mg"* over a Hubble time. 

It should be noted that this explanation of the first Sx problem differs somewhat from those previously proposed (see 
§2.2). In particular, Blakeslee (1997) and Blakeslee et al. (1997) showed a correlation between the specific frequencies 
on 40-kpc scales in BCGs, and the X-ray luminosities over much larger volumes {Rgc < 500 kpc) in their host clusters. 
These authors, along with West et al. (1995) and Harris et al. (1998), therefore interpret the phenomenon of high 
specific frequency as one associated in some way with entire clusters of galaxies. The analysis presented here, however, 
has attempted to address the problem more locally, by referring to the gas and stellar masses on similar spatial scales 
(typically rgc <; 100 kpc) to estimate G (eq. [^) in individual systems; nowhere has any appeal been made to more 
global properties of galaxy clusters. This general approach is motivated, of course, by the observed situation in M87 
specifically (§3). The correlation that does exist between BCG Sx and cluster- wide Lx (and other such relations; see 
the references cited) may simply stem from the fact that the gas around BCGs in rich clusters blends smoothly into 
the intracluster medium as a whole. Note that a relation of some kind between "local" Lx and Atot is also implied by 
Santiago & Djorgovski (1993), although they attach no specific interpretation to it. 

Interestingly, the bold line in Fig. H departs significantly from linearity (i.e., from a simple power-law scaling between 
A/tot and -Z^ygai) for Mgas/Mgtars ^ 0.1-0.2, corresponding to Ly > 3-5 x lO^^L© or My < —21.5. This magnitude is 
similar to that which appears to divide ellipticals into two broad classes in terms of their kinematics, isophote shapes, 
and light distributions. (Lower-luminosity ellipticals generally show rapid rotation, fiattened and disky isophotes, and 
steep cusps in their central densities, while brighter systems are slowly rotating, with round and boxy isophotes and 

^''The term BCG is used loosely here, to refer to the centrally dominant galaxies in groups as well as in rich clusters, even when these are 
not strictly ranked first by luminosity. 
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Fig. 8. — Observed total GCS populations, A/tot, vs. total galaxy luminosity, Lv,gal! for a sample of 65 giant ellipticals and BCGs (filled 
and open circles; from Harris et al. 1998, and references therein). Dotted lines have A/tot oc Lyg^^i, or constant Sm = 15, 5, and 1.5. The 
systematic departure of BCGs from Sjv ~ 5, at I/v^gai > ~ 5 X 10^*^ Lq, is an illustration of the first specific frequency problem. The bold, 
solid curve is the prediction of an assumed universal cluster formation efficiency e^i = 0.0026; the thin, dashed line is what that prediction 
would be if all ellipticals were completely devoid of gas. The two most extreme outliers are NGC 3557 and NGC 5018, each of which is 
relatively isolated and has Ly ~ 1O^"'^L0 (My ~ —22.5) but Sm < 1- 



shallower central density profiles; see, e.g., Davies et al. 1983; Bender et al. 1989; Nieto, Bender, & Surma 1991; Tremblay 
& Merritt 1996; Gebhardt et al. 1996; Faber et al. 1997.) This bears directly on the recent discussion of Kissler-Patig 
(1997; see also van den Bergh 1998), who used observed trends in GCS specific frequencies and radial distributions as 
further evidence for the empirical "dichotomy" of ellipticals. In particular, Kissler-Patig noted that the cluster systems 
in early- type galaxies with My <; —21.5 tend to have lower specific frequencies and steeper A'ci profiles than those in 
brighter galaxies. Although he went on to argue for a discontinuous change in the cluster formation efficiency for galaxies 
with Ly below and above « 5 x 10^^ Lq, Fig. |^ now shows that all of the data are well represented by a continuously 
increasing Sn that derives from a constant ed and a smoothly increasing gas mass fraction, G; the nonlinearity intrinsic 
to the resulting Atot^^v.gai relation just happens to roughly mimic a break of sorts around My — —21.5. 

This interpretation can also account for another correlation suggested by Kissler-Patig (1997), namely, that Sn appears 
to be lower in (faint) ellipticals whose isophotes show disk-like deviations from pure ellipses, and higher in those (brighter) 
galaxies with more box-like shapes. These isophotal perturbations are themselves correlated with the X-ray luminosities 
of ellipticals (Bender et al. 1989), in the sense required qualitatively to explain this observation of GCSs again in terms 
of a constant ed and a systematically increasing gas-to-stellar mass ratio. This raises the interesting question of which, 
of gas content or isophote shape, is the more fundamental property of elliptical galaxies; the present discussion, at least, 
would appear to favor the first alternative. 

Finally, although the slopes of GCS radial distributions cannot be addressed directly by Fig. ||, it seems natural that, 
given a fixed ed, ellipticals with large gas contents (and thus high global Sn) should generally have shallower GCS density 
distributions, relative to their stellar profiles, at large galactocentric radii (simply because Pgas/Pstars usually increases 
with Tgc in early-type galaxies). As was suggested in §2.2, then, the first and second Sn problems should be closely 
related, with G = Mgas/Mgtars the controlling factor in both. Again, this is consistent with the observed trend (noted by 
Forbes et al. 1997; Kissler-Patig 1997; and van den Bergh 1998) towards weaker dependences of A^ci on i?gc in brighter 
and higher-S'jv galaxies. This proposal should, of course, be checked on a case-by-case basis, as it was for M87 in §3.2 
above, in as many systems as possible. 

4.2. Dwarf Ellipticals 

If the bold line in Fig. || is extrapolated into the regime of dwarf ellipticals, Ly.gid < 2 x lO^L©, it clearly implies that 
global Sn < 2 should be observed in these faint galaxies. This is completely at odds with the data of Durrell et al. (1996) 
and Miller et al. (1998), who find an average Sn « 5 for a large sample of dE's in the Local Group, in other groups 
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Fig. 9. — A/tot vs. iv.gal for early-type systems, now including 32 dwarf ellipticals with Lv.gal ^ 2 X W^Lq (solid squares; Durrell et 
al. 1996 and Miller et al. 1998). Filled and open circles represent the same giant ellipticals and BCGs plotted in Fi|ft. pi The dotted lines of 
constant 5jv = 15, 5, and 1.5, and the heavy solid curve with e^l = 0.0026 are also the same as those shown in Fig. M. The bold, dashed line 
is the prediction for the same constant e^l = -'^g?s'/^^gas' = 0.0026, after correcting for significant gas loss in supernova-driven winds from 
dwarfs fainter than My ~ —18.5. The large open square is plotted at the observed Ly and A/tot of the Milky Way spheroid; see §4.3 of the 
text. 



nearby, and in the Virgo and Fornax Clusters. In addition, Miller et al. claim evidence for an increasing Sn towards 
fainter i\/,gai- 

This can be explained by first recalling that dwarf ellipticals do not fall on the fundamental plane of bright galaxies; 
thus, the scalings (^5|) and (|2^) do not apply to them. Also, these galaxies are bluer than their giant counterparts, so 
that the conversion to y-band mass-to-light ratios from measured T^j, as described above, is again inapplicable; and 
in any case, even the cores of dwarfs appear to be dominated by dark matter, so that their dynamical Ty are not 
necessarily representative of the stellar values. Most importantly — and likely at the heart of these other points — dE's 
present sufficiently shallow potential wells that they should have suffered significant amounts of mass loss during supernova- 
driven winds in early bouts of star formation. Thus, even though these galaxies are gas-poor at the present, the appropriate 
observable Cci = Mgcs/Mgtars is not an accurate estimator of the true cluster formation efficiency ed = M^'^^ / (see 
Durrell et al. 1996 for the first suggestion that this is an important consideration in understanding dwarf GCSs). Even 
still assuming that the GCSs now have Mgcs — M™]}", allowance has to be made for an additional factor Mstars/Afi^s* < ^■ 
The relation analogous to equation ( p^ ) in this case is therefore 

4 Mgcs / ^^gas* ( f Lv,gal \ ( Ty \ f (m)cl 



- 1-7 X 10 ^ j [j^j^) [-^J^ j (,2.4x10^mJ ' 

where Ty ~ 2 should be a good representation of the stellar populations in systems fainter than about 2 x 10^ Lq. 

The specific model of Dekel & Silk (1986), which considers supernova-driven winds in dwarf galaxies that are dynamically 
dominated by dark matter halos, can be used to evaluate the ratio of initial gas to present stellar mass. That study finds 
that the winds from systems with velocity dispersions less than about 100 km s~^, corresponding to My ^ —18.5 at 
the current epoch, should have expelled whatever gas remained after a single burst of star formation. Moreover, as one 
moves towards shallower potential wells and fainter Ly.gai, the fractional mass lost, relative to the initial total, is naturally 
expected to increase. The predicted scaling, in the limit of very large mass loss and for a power spectrum similar to that 
in a standard cold dark matter cosmology, is roughly Afg"g*/Mstars oc iga"'"'- Taking the crude step of extrapolating this 
all the way to the critical velocity dispersion (i.e., to iy,gai — 2 x lO^L©), and assuming that any gas driven from brighter 
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and more massive galaxies was a negligible fraction of the total, the implication is that 



gas 



2 X 109 Lq 



^0 



(28) 



Substituting this in equation (|27|), and assuming a constant stellar Ty — 2, yields A/tot oc iygai and 5jv oc iygai for 
dwarf ellipticals. Although the picture of dE's as single-burst populations is now known to be incomplete (see, e.g., Mateo 
1998), this result is nevertheless consistent with the observations of Miller et al. (1998). 

Figure |9| uses the data from Durrell et al. (1996) and Miller et al. (1998) to place 32 nucleated and non- nucleated dwarf 
ellipticals into the Lv,gai — A/tot plane with the brighter galaxies from Fig. ||. (Of these, 16 are taken from Durrell et al., 
and 16 from Miller et al. Only GCSs detected at the > l-tr level are used.) The bold dashed line in the lower left of this 
plot comes from equations (p7| ) and (p8|), assuming a constant ed = Mgcs/Mg"g* = 0.0026. The heavy solid line is the 
same as that in Fig. ^ — i.e., it also has ed = 0.0026, but applies to fundamental-plane ellipticals — and the light, dotted 
lines correspond, as before, to constant Sn = 15, 5, and 1.5. Once again, all the data appear consistent with a universal 
globular cluster formation efficiency. It will be important to properly confirm this with detailed models for the dE's in 
particular, and to flesh it out with observations of systems at iv,gai '--^ (1 — 5) x lO^L© in general. 

4.3. Globular Clusters in the Milky Way 

One object that falls just inside this observational "gap" in luminosity is the stellar spheroid — disk excluded — of the 
Milky Way. The total F-band luminosity here has been estimated by de Vaucouleurs & Pence (1978) to be Ly — 
5.1 X lO^L© [after applying a color of {B — V) = 0.65 to their Mb — —18.8]. At the same time, the total number of 
globular clusters that are genuine halo objects, rather than part of the (thick) disk, is Atot — 110-120. (Halo globulars 
are identified on the basis of their metallicity, [Fe/H] < —0.80; refer back to Fig. ^ above.) This is the point plotted as 
the large open square in Fig. ||. 

In fact, de Vaucouleurs & Pence (1978) model the Galactic spheroid with an R^^^ density law with effective radius 
i?cff = 2.67 kpc. Thus, in addition to predicting the absolute magnitude that would be seen by an external observer, it is 
also possible to parametrize the density profile of the entire stellar halo. To do this in three dimensions, Young's (1976) 
deprojection of the i?^/^ law may be used: 



Pstars (^gc) 
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(29) 



g and i?off directly from the work of de Vaucouleurs 



which is vahd for rgc ^ 0.2i?off- Substitutions can be made for L^* 
& Pence. One way to estimate the B-band mass-to-light ratio of the Galactic Population II is to adopt the same T-Lb 
scaling used in §4.1, but to normalize the relation by adopting solar colors rather than the redder ones that apply in 
ellipticals, van der Marel's (1991) mean = 3.32/i5o then corresponds to Tb — 4.65M0Lq^ at Mb — —21.49 (for 
Hq = 70 km s^^ Mpc^^ assumed), and cx i^"^. The implied TB(halo) ~ 2.2MqLq^ for Mb — —18.8 is consistent 
with the mean Ty = 2 for globular clusters, and with the Ty = 1.75 adopted for the spheroid by Bahcall & Soneira 
(1980) in their model of the Galaxy. Equation (|29| ) is therefore evaluated here with Tb = "^MqLq^ . 
The dashed line in the bottom panel of Fig. |2| in §3.1 above shows this density profile, scaled down by a factor of 



;mw 



= 0.0027 ±0.0004 



(30) 



This value is obtained as the mean of the ratio Pd/Pstars over all galactocentric radii 2 < rgc < 40 kpc, and it accounts for 
the position of the Milky Way halo in Fig. ||. Such a basic correspondence between our Galaxy and early Hubble types 
was previously suggested, although developed from a somewhat different argument, by de Vaucouleurs (1993). 

It is also evident from Fig. ^ that there is no second Sn problem in the Milky Way: pci oc pstars obtains throughout 
the entire halo at rgc ^ 3 kpc — so again, beyond roughly an effective radius, within which the original GCS might have 
suffered significant dynamical depletion. This point has long been appreciated, although it is usually discussed in the 
context of the fact that the same rough rule p oc r^J^-^ is inferred both for the spheroid, from counts of RR Lyrae and 
horizontal-branch stars (e.g., Preston, Schectman, & Beers 1991; Kinnian, Suntzeff, & Kraft 1994), and for the GCS, from 
simple power-law fits to its density profile (Harris 1976; Zinn 1985; Djorgovski & Meylan 1994). This power law is shown 
as the dotted line in the bottom panel of Fig. |^. 

4.4. Scatter in the Ly-Atot Plane 

As a final test of the general viability of a single globular cluster formation efficiency, Fig. [l^ plots the ratio of observed 
to predicted GCS Atot as a function of host galaxy luminosity for all of the 98 systems discussed thus far. The predicted 
Atot for fundamental-plane ellipticals and the Milky Way spheroid (filled circles and open square in the Figure) are given 
by equation (p^, with the gas mass fraction G set identically to zero; for BCGs (open circles), by the same relation 
but with G as in equation (p6[); and for dE's (filled squares), by equation (27). In all cases, it has been assumed that 
ed = 0.0026. The median of the resulting ratios 7^ = Af^^^" /M^^^'^ is 1.04, and the geometric mean is 0.98. Close to 80% 
(77/98) of these galaxies have ^ < 7^ < 2. 
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Fig. 10. — Ratio of observed Mtnf to the predicted values for all GCSs represented in Fig. |9| using e^l = 0.0026 in equation ( [24| ) or 
Mgca/M^n^* = 0.0026 in equation (|^. Point types have the same meaning as in Fig. |. The ratios for NGC 3557 and NGC 5018 are 0.07 
and 0.17, and do not appear on the scale of this plot. 



As was also mentioned in §4.1, but can perhaps be seen more clearly here than in Fig. ^ or Fig. ^ the deviations 
of observed A/tot from the uniform-ed prediction show no significant dependence on either galaxy luminosity or Hubble 
type. There are possibly two exceptions to this claim: (1) the GCS populations in the brightest BCGs, -/jy,gai ^ 10^^ Lq, 
appear to fall below the values predicted with ed = 0.0026, by perhaps ^ 50% on average; and (2) the two faintest 
dwarfs, which happen to be the Local Group sphcroidals Fornax and Sagittarius, both have A/tot greater than predicted. 
These rough impressions may well not remain, in the mean, after future additions of GCS data from other very high- and 
low-luminosity galaxies; but even if they do, they need not imply the existence of fundamental variations in ed. Rather, 
there is likely enough uncertainty in the adopted relation between gas mass fraction and Ly.gai (eq. [ p6| ) to account for 
the first item; and the second could reflect either a deficiency in the very simple correction (§4.2) for gas blow-out from 
very faint dwarfs, or the neglect of possibly extensive post-formation dynamical evolution in such low-mass galaxies (such 
as tidal stripping of their field stars; recall, for example, that Sagittarius is in the process of being accreted by the Milky 
Way). 

Even if these issues are ignored completely, a least-squares fit to all of the data shown in Fig. reveals an exceedingly 
weak dependence oi Af°^^ / Af^^^'^ on galaxy luminosity: TZ cx gai**°'°'^- Fitting to only the fundamental-plane ellipticals 
and BCGs yields a steeper TZ cx ij^^fi or so, but this result appears to be driven by just three gE's with anomalously 
low or high GCS populations (cf. Fig. |): NGC 3557, with Ly = 9.4 x 10^° Lq and Sn = 0.4 ± 0.3; NGC 5018, which 
has Lv = 7.8 x 10^° Lq and Sn = 0.9 ± 0.3; and NGC 4278, with Ly = 7.1 x 10^ Lq but Sn = 12.3 ± 1.4. If these 
systems are excluded, then 7^(gE + BCG) cx Ly°^^^'°^ is indicated instead. Similarly, the slopes of TZ vs. Ly.gai do not 
differ significantly from for any of the individual Hubble-type samples represented in Fig. |l^. 

Thus, current observations of total GCS populations do show some scatter, and it will be important to understand — 
likely through detailed analyses, along the lines of §3, of more systems individually — whether or not this reflects an 
intrinsic scatter in the underlying globular cluster formation efficiency (see also §4.1). However, there is no clear evidence 
at present for any systematic departures from a mean global cfhciency of ed « 0.0026. 

4.5. Open Clusters 

As was discussed in §2.1, the formation of stellar clusters continues to be an important element of star formation 
in general, and the existence of an apparently universal formation efficiency for globular clusters must reflect a generic 
piece of this process. However, the picture may still be incomplete: although the discussion to this point has included 
clusters that formed in a wide variety of environments — at galactocentric distances rgc ~ 10-100 kpc, in galaxies spanning 
four orders of magnitude in total luminosity and inhabiting both poor groups and rich clusters — all of these formed at 
essentially a single epoch, and under physical conditions that were possibly much more extreme than those in, say, the disk 



Mclaughlin 



27 



of the Milky Way today. It is important, therefore, to compare the protogalactic ed — 0.0026 ± 0.0005 to the frequency 
with which bound star clusters form in a typical quiescent galaxy at the current epoch. Disk clusters in the Milky Way 
present one opportunity to do this. 

The rate of open cluster formation in the solar neighborhood has been estimated by Elmegreen & Clemens (1985) as 
Sci — (2.5 ± 1) X 10~^ clusters kpc~^ yr~^- If the masses of the clusters range between 100 M0 and 5OOOAf0 according 
to dM /dm cx m~^-^ (from the luminosity spectrum d/V /dL of van den Bergh & Lafontaine 1984), then a mean mass of 
~ 700 M0 is indicated, and the formation rate becomes Sc\ — (1.8 ± 0.7) x kpc^'^ yr^^. This is to be compared 

with the total star formation rate in the solar neighborhood, which includes those stars appearing in loose groups and 
unbound associations (or even in relative isolation) as well as in bound clusters. This rate is 5tot — 4.5 x 10"'^ Mq kpc~^ 
yr-i (e.g., McKee 1989). 

Since the mean stellar mass is no different in clusters than in the field, the ratio of these two rates can be taken directly, 
to assess the local number fraction of stars born in clusters: Sc\/Stot ^ 0.04. This is consistent with the classic estimate 
(e.g., Roberts 1957) of about 10%. However, the globular cluster EcI has been obtained here as a mass ratio, since this 
is the observable quantity that can be related directly to its formation value (see eq. [^). To make a proper comparison 
with open cluster formation, it is therefore necessary to refer the mass going into newborn clusters (per unit area, per 
unit time) to the total mass of gas available for any kind of star formation. This latter figure can be derived from iStot 
using the mean star formation efficiency, (SFE) = (Mgtars/Afgas), averaged over many of the dense, star-forming (not 
necessarily cluster-forming) cores in giant molecular clouds. This number is only poorly known, but statistical arguments 
(e.g., Elmegreen 1983) suggest that it lies at the level of (SFE) - 5%-10%. Thus, 

eci(open) ~ s^^^ff^^^^ - 0-002 - 0.004 , (31) 

which is in reasonably good agreement with the result for globular cluster systems. The numerical value in equation 
( ]3l| ) clearly should not be taken as definitive, given the very rough, order-of-magnitude nature of the argument that has 
produced it. That said, however, it is just the order-of-magnitude agreement between equations ( ^l| ) and ( p3| ) which 
suggests that the basic mechanism of cluster formation may not have changed appreciably from protogalactic times to 
the present. 

Another obvious place to check the current level of ed is in starbursts and merging galaxies, which are forming massive 
young clusters in substantial numbers (e.g., Whitmore & Schweizer 1995; Meurer et al. 1995). Complete inventories of the 
star-forming gas mass in individual systems, and of the total mass in those of their young clusters that can be shown to 
be gravitationally bound (which to date has been done rigorously for only two objects: Ho & Filippenko 1996a, 1996b), 
would be of considerable interest. The results should be of particular relevance to the oft-made claim, which is now at the 
center of some debate (Sternberg 1998; Brodie et al. 1998), that these super star clusters are the modern-day equivalents 
of young globulars. If it turns out, for example, that strongly bound clusters in starbursts are forming in numbers very 
much in excess of those expected on the basis of eci — 0.25% (as, e.g., the UV observations of Meurer et al. 1995 appear 
to suggest), then these spectacular events might not be representative of a ubiquitous phase in the formation of average 
elliptical galaxies. 

5. DISCUSSION 

The fact that the ed values inferred here are everywhere so similar — over large ranges in radius inside M87, M49, and 
NGC 1399; from dwarf ellipticals to brightest cluster galaxies; and from massive, old globular clusters to smaller, much 
younger open clusters — should serve as a strong constraint on theories of star formation. The simplest interpretation is 
that the probability of attaining a cumulative star formation efficiency of SFE ^ 20%-50% (see §2.1) in any dense clump 
of gas more massive than ~ 10^-10'^ Mq (so from open clusters on up) can depend only weakly, on average, on local 
environment (i.e., the density and pressure of ambient gas) or on factors such as global background density (meaning, e.g., 
in the protogalactic context, a l-cr vs. a 3-i7 density fluctuation; or a poor group of galaxies vs. a rich cluster). Moreover, 
implicit in plots like Figs. |^ and|^ is the fact that cluster formation efficiencies are also rather insensitive to metallicity: the 
mean abundances of the globulars in the galaxies represented there range over two orders of magnitude, from [Fe/H] ~ —2 
to [Fe/H] ~ (see Durrell et al. 1996; Forbes et al. 1997). Over a remarkably broad spectrum of physical conditions, 
therefore, it appears that similar fractions of star-forming gas — always about 0.25% by mass — are able to produce bound 
stellar clusters. The implicitly mass-averaged nature of the observational estimates of ed has to be kept in mind, but 
would seem to be of little concern in most cases. It is therefore reasonable to place the efficiency of cluster formation 
alongside the globular cluster mass function dN /dm (McLaughlin & Pudritz 1996) and the mean globular cluster mass 
(m)d (e.g.. Fig. ^ above) as robust and nearly universal physical quantities that must be derived from the most general 
aspects of the star formation process. 

As a corollary to this, one implication of Figs. ^ and ^ above is that the efficiency of unclustered star formation could 
not have been universal. That is, in both dwarf elliptical and brightest cluster galaxies, where the data are consistent 
with ed = .Wgcs*/-^gas' having been constant in the mean, globular clusters apparently did form in precisely the numbers 
expected for the total amount of gas initially available to these systems. But then, the higher than average specific 
frequencies there are inferred to reflect larger ratios of M™!,*/Mstars, so that lower than average amounts of this gas were 
converted into field stars. The data for dwarfs agree with a simple scenario in which most of the unused gas supply is 
driven entirely out of the galaxies by a strong wind following one major burst of star and cluster formation. If this is 
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basically correct, then the globulars in dE's quite possibly were able to form completely before the onset of such a wind, 
and they could conceivably have been instrumental in driving it; at the very least, if the protoclusters were still largely 
gaseous, they must have been already sufficiently dense and well defined that they could survive such a catastrophic 
event essentially intact. At the same time, field star formation must have been truncated by the wind, before reaching 
what would have been a "normal" efficiency of Mstais/M^^^ = 1 — ed- Thus, the massive and dense clumps of gas which 
ultimately formed bound star clusters had to have collapsed more rapidly than those which produced unbound groups and 
associations. This conclusion may hold quite generally — one possible interpretation being that only those gas clouds in 
the highest-density tail of some distribution, in any setting, are able to achieve a high cumulative SFE — and it may be 
related to the fact that the globular clusters in dwarfs are generally somehat bluer, and presumably more metal-poor, 
than the bulk of the field stars there (Durrell et al. 1996; Miller et al. 1998) 

A similar feedback argument could also apply to BCGs, where it may be that especially strong early bursts of star and 
cluster formation led to the premature virialization of large amounts of protogalactic gas, and perhaps drove slow, partial 
galactic winds, most effectively in the low-density environs at large galactocentric radii (cf. Harris et al. 1998). The unused 
gas would have to remain hot to the present day, and more or less in the vicinity of the parent galaxy, in order to enter 
the observed Lx-Lb relation as inferred in §4.1. This requirement is consistent with the fact that BCGs are situated at 
or near the centers of their galaxy clusters, and thus at the bottoms of very deep potential wells. In this scenario, again, 
globulars would have had to form, in just the numbers expected of them, somewhat before field stars were able to do the 
same; and this truncation of unclustered star formation would have had to be more severe at large galactocentric radii 
(where the feedback is strongest). This could then account explicitly for the connection between the first and second 
specific frequency problems, i.e., between high global Sn and local Sd/Sgtars ratios that increase with i?gc. It also gives 
a more specific context to the general claim (Harris 1986) that the shallowness of CCS radial distributions relative to the 
stellar light profiles in some ellipticals resulted from the globulars having formed slightly in advance of the stars there. The 
main demand on this picture for BCG formation is that it must be able to explain why Mgas/Mgtars increases with galaxy 
luminosity, i.e., why more massive BCGs apparently had larger fractions of their initial gas mass virialized early on. As 
Harris et al. (1998) discuss in some detail, this issue is likely related to the fact that brighter BCGs are generally found 
in more massive galaxy clusters, which may have presented higher-density and more turbulent environments that led to 
systematically higher star formation rates, and more violent feedback. Although speculative, these ideas are suggested 
directly by the CCS data, and are consistent with all those available. Since they are also obviously related to questions on 
the star-formation histories of dE's, on the origin of the Lx-Lb correlation in the brightest ellipticals, and on the nature 
of the intracluster medium, they should be addressed with considerably more rigor than has been applied here. 

It is particularly striking that the cluster formation efficiency derived in §4.3 for the Milky Way halo is essentially 
indistinguishable from that which applied in M87, M49, and NGC 1399, and which appears to have held quite generally 
in ellipticals of most any description. This calls into serious question one of the primary motivations for the model of 
Ashman & Zepf (1992) and Zepf & Ashman (1993), who posit that ellipticals are formed primarily by the mergers of 
gas-rich spirals, and invoke preferential cluster formation during these events to account for what they claim is a typical 
factor of two difference in the observed ratios Mgcs/Mgtars of late- and early-type systems. Although it is true that 
spirals have specific frequencies that are lower than those in bright E galaxies, by factors of perhaps 2-3 on average (as 
is exemplified by the Milky Way in Fig. ^ above; see also Harris 1991), it is now apparent that this is not necessarily a 
reflection of similar discrepancies in their more fundamental Cci- If the Milky Way is typical, it suggests that — even though 
stellar clusters certainly can form in merging galaxies — mergers are not required to explain the relative GCS populations 
of spirals and ellipticals in general. Conversely, whenever an elliptical does form by the major merger of two spirals, this 
is allowed to include cluster formation only at the standard efficiency of ed — 0.0025. 

Ashman & Zepf (1992) and Zepf & Ashman (1993) also attribute the second specific frequency problem to the formation 
of ellipticals by mergers of spirals with pre-existing GCSs. However, it is not clear that this is necessary either, if the 
extended spatial distribution of globulars in the brightest ellipticals is generally related to their association with hot gas 
at rather large galactocentric radii. In fact, a basic expectation of the original Ashman-Zepf scenario is that the second 
Sn problem should be less pronounced in galaxies with higher global Sn (which they hypothesize to have suffered greater 
numbers of major mergers). The sense of this predicted trend is opposite to the one implied here — with empirical support 
directly from the evidence in M87, and indirectly from the correlations of Forbes et al. (1997) and Kissler-Patig (1997) — in 
which a constant Cci combines with the larger Afgas/A/gtars in higher-S'Ar ellipticals, to produce greater contrasts in their 
projected Sci and Sstars profiles. Although radial-profile and Sn data by themselves may not rule out the merger model 
for GCSs altogether, a revision to take account of the hot gas in large ellipticals is essential if it is to be at all viable. (See 
also Forbes et al. 1997 and Kisslcr-Patig et al. 1998 for critical discussions of this and other aspects of the basic idea.) 

Finally, the GCS density profiles derived specifically for M87, M49, and NGC 1399 in §3 also are relevant to issues of the 
generic spatial structure of GCSs, and to the specific matter of the possible existence of intergalactic globulars in clusters 
of galaxies. To facilitate a discussion of these points, and for reference. Table || presents the results of least-squares fits of 
a few different parametric functions to the spatial distributions of the M87, M49, and NGC 1399 GCSs, as they are given 
in Tables |2[ H, and ^ above. The functions include (from left to right in Table ||) a power law with a constant-density 
core region; the cusped density profile suggested by Navarro, Frenk, & White (1996, 1997) as appropriate for cold dark 
matter halos; and the law popularized by Hernquist (1990) and applied widely to the light profiles of elliptical galaxies. 
In all cases, fits were performed by minimizing the residuals of the functions' projected density profiles relative to the A^ci 
data, which have smaller observational uncertainties than the volume densities do. The best-fit parameters thus obtained 
were then used to overlay the original, three-dimensional functions on the deprojected and verify that they provided 
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good descriptions of those as well. One consequence of this procedure is that no fit is given for the Hernquist (1990) 
profile, rici cx; (r/a)^^(l + r/a)"^, to the M87 GCS. Although it is possible to fit this function directly to the n^i data, 
its projection affords a poor description of the full A'ci distribution. This almost certainly reflects the fact that the GCS 
profile in M87 appears about to become significantly shallower, at least temporarily, beyond the last radius (rgc ~ 100 
kpc) for which cluster counts exist (because Fig. || suggests essentially that pd cx pgas there). Such a change in slope is 
obviously not allowed by Hernquist's density profile. Thus, the failure of this model in this case makes the point that any 
of the density fits in Table |^ should be taken at face value only over the ranges of rgc that have been directly observed; 
extrapolations to much larger radii could potentially introduce appreciable systematic errors. 

Table |]has three interesting implications: (1) None of the available data require the existence of truly constant-density 
"cores" in GCSs. This was also mentioned by McLaughlin (1995) and it is, of course, plainly seen in Figs. ^ through p[ 
Evidence is emerging for a similar situation in elliptical galaxies generally (Grane et al. 1993; Gebhardt et al. 1996), and 
apparently also in galaxy clusters (Merritt & Tremblay 1994; Carlberg et al. 1997; McLaughlin 1999). (2) It has been 
suggested, particularly in the contexts of some dynamical analyses (e.g., Huchra & Brodie 1987; Weil, Bland-Hawthorn, 
& Malin 1997) that the M87 GGS directly reflects the distribution of dark matter in and around that galaxy. However, 
this is not true in any simple sense; neither of the rid fits given in Table || provides an adequate description of /3dark(?'gc) 
in M87. Instead, it is found there that pdark oc pgas (see McLaughlin 1999), in which case Fig. ^ shows that the GCS 
becomes a reliable tracer of the dark matter halo alone only on fairly large spatial scales (where the number densities 
become quite low). (3) The analytical fits to A''ci and offer a convenient way to extrapolate the observed densities to 
(slightly) smaller and larger radii than those observed, so as to compute the total GGS populations of these three galaxies. 
Thus, in M87, it is found that A/tot = 13 600±500 clusters (over all magnitudes) are projected to distances Rgc < 25' ~ 110 
kpc from the galaxy's center; in M49, there are Atot = 6 850 ± 550 globulars within the same projected Rgc, and in NGC 
1399, Atot = 4 700 ± 400 clusters have Rgc < 14'5 ~ 70 kpc. These total "projected populations" are in good agreement 
with those estimated in many other studies (see the references cited in §3). However, it is now possible to determine the 
three-dimensional radius of the sphere which actually contains all of these clusters in each system. From either of the M87 
fits in Table ||, it is found that 13 600 clusters are contained within rgc — 150-160 kpc; in M49, 6 850 globulars come from 
the sphere rgc < 140-150 kpc; and the 4 700 globulars in NGG 1399 all have true galactocentric positions rgc < 90-95 
kpc. To put this another way, all three GCSs are sufficiently centrally concentrated that, in each case, some 85%-90% of 
the globular clusters projected onto the galaxy — so again, to within i?gc < 110 kpc for M87 and M49, and Rgc < 70 kpc 
for NGC 1399 — really are located at those galactocentric distances in three dimensions. 

This last point specifically can clarify and constrain the suggestion (White 1987; West et al. 1995; Cote et al. 1998) that 
the high specific frequencies in BCGs like M87 are due to the "contamination" of their GCSs, in projection, by a significant 
number of globular clusters that are associated with a galaxy cluster as a whole, rather than bound to the central galaxy 
itself. Harris et al. (1998) tried to test this idea in some detail for M87 in particular, working from the assumption 
(following White 1987) that any population of intergalactic globular clusters (IGCs) should generally be associated with 
the diffuse stellar light of the cD envelopes that surround many BCGs (i.e., these structures were considered to comprise 
primarily intracluster stars). This hypothesis leads to inconsistencies that prompted Harris et al. to reject IGCs as the 
main cause of the first specific frequency problem. While their arguments still stand within the framework that they 
constructed, it is now seen that the high Sn of M87 certainly, and likely those of other BCGs as well, can be traced to 
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"'^For radii in arcminutes and densities in number per arcmin^. 

^Normalizations no apply to y < 23.9; divide by 0.5572 to obtain population over all magnitudes. 
•^Normalizations no apply to V < 23.3; divide by 0.3646 to obtain population over all magnitudes. 
""^Normalizations no apply to V < 24.0; divide by 0.5332 to obtain population over all magnitudes. 
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subsets of their GCSs that are associated with hot gas, and not necessarily with stars at aU. This does appear to favor the 
basic notion of IGCs, since the gas around M87 extends to very large spatial scales, and traces out the potential well of 
the Virgo Cluster as a whole (McLaughlin 1999). If a constant eci holds all the way out, then, a cluster- wide population of 
globulars could be expected. However, since all of the 13 600 globular clusters which give M87 its high specific frequency 
in projection are physically associated with a volume of radius ~150 kpc around the galaxy, it is not obvious that any of 
them are truly "intergalactic" objects at the current epoch. Projection effects arc not solely responsible for the high S* Az- 
of M87, and although intergalactic globular clusters may well exist, it is not clear that any have yet been observed. 

Of course, since M87 is at the dynamical center of Virgo, it can be difficult to distinguish meaningfully between the 
galaxy and the cluster on 100-kpc scales. In the context of IGCs, the question evidently comes down to one of origin: 
did the globulars now seen right around M87 actually form there, or were many of them somehow brought in from 
much further than 150 kpc away? Any self-consistent answer to this must be able to account for the observed fact that 
Pel = 0.003 (pgas + Pstars) at galactoccntric radii r^c ~ 10-100 kpc in M87; that is, both the very existence and the 
normalization of this proportionality have to be explained. Although essentially any scenario which is able to do this 
would be allowed by the current data, none could be preferred over what remains the simplest option: that the bulk of the 
M87 GCS formed in situ, with a thoroughly standard and spatially constant EcI- Ultimately, only detailed evolutionary 
models, tailored specifically to the properties of M87 and Virgo and considering both stars and gas, can really decide 
what fraction of the galaxy's GCS might have come originally from the cluster at large. 

6. SUMMARY 

Stars form mainly in groups, some of which emerge from their natal clouds of gas as gravitationally bound clusters. 
An understanding of cluster formation will therefore be integral to any full theory of star formation, and an important 
element of this is simply the frequency with which it occurs. That is, what is the likelihood that a star-forming cloud of 
gas whether this be a dense clump in a Galactic giant molecular cloud, or a larger one in a protogalactic fragment will 
achieve the high cumulative star formation efficiency (viz. SEE ^ 20%-50%) that allows its stars to remain bound as a 
cluster after they have cleared away the gas? This paper has used observations of the globular cluster systems in galaxy 
halos to address this question, and to empirically evaluate the efficiency of cluster formation: by mass, ed — 0.25%. This 
is, by all appearances, nearly a universal number, and it should therefore serve as a strong constraint on theories of star 
and cluster formation in any context. 

To arrive at this result, it was first shown (§2.1) that the total (global) masses of GCSs are quite robust, even over 
a Hubble time, against the dynamical processes that work to destroy low-mass globular clusters in the central regions 
of galaxies. With an observed Afgcs thus useful as a measure of the initial quantity, the specific frequencies of GCSs, 
^jv oc A/tot /-^v, gal c< Mgcs/Mgtars, Serve as crude first estimates of the basic formation efficiency ed = -^gcs'Z-^gas*- To 
the limited extent that 5jv ~ 5 is an adequate description of early-type galaxies, the implication is that ed ~ 2 x 10~^. 
The limiting factor in this simplest treatment of GCS data is the assumption that a galaxy's current stellar mass always 
suffices as an estimate of its total initial gas mass; if taken too far, this quickly leads to some perplexing conclusions. 

As was reviewed in §2.2, systematic variations in Sn exist along the entire sequence of early-type galaxies, from dwarf 
ellipticals (where A/tot oc Ly^^^^ or so) through normal giants (A/tot oc g^i) and on to the brightest ellipticals, including 
many of the central galaxies in groups and rich clusters (Atot oc -'^y gai)- If specific frequency were indeed a good estimator 
of ed in general, the implication would seem to be that this also varied, for unknown reasons and in a non-monotonic 
fashion, as a function of galaxy luminosity. Closely related to this is the fact that the GCSs of some, but not all, large 
ellipticals are spatially more extended than their stellar halos. That is, local 5jv values increase with projected Rgc in 
some systems, which naively would suggest that — again, for unknown reasons, and counter to any intuition — the efficiency 
of cluster formation sometimes increased with galactoccntric radius. 

An important clue to the true nature of both these "specific frequency problems" which are most likely two aspects 
of a single phenomenon — lies in the observation that the global S'jv in central cluster galaxies increases with the X-ray 
luminosity of their parent clusters. Some of the globulars in these systems ought then to be associated with the hot gas 
there; and is approximated better by the present-day (Mgas +Mstars) than by Mgtars alone. Sn then systematically 

overestimates the true cluster formation efficiency ed, which instead is given empirically by Cci = -^gcs/(-^gas + -^stars)- 
In addition, since a galaxy's hot gas generally follows a shallower density profile than its stellar halo does, the "extra" 
clusters in very gas-rich (and high-S'Tv) ellipticals might well be expected to also show a radial distribution, pd or Ed, 
that is distended relative to Pstars or 

These ideas, which extend similar suggestions in the recent literature, have been tested in detail for each of M87, M49, 
and NGC 1399. Data were collected from the literature on the GCSs, surface brightness profiles, and X-ray gas densities 
of these galaxies, covering spatial scales rgc ~ 1-100 kpc. In §3.1, new comparisons of the projected densities Ed and 
Estars in M49 and NGC 1399 (which are gas-poor) showed that, contrary to previous claims, neither galaxy suffers from 
either of the two specific frequency problems: The global Sn of M49 has a "standard" value of ~5, and the local ratio 
Sd/Egtars is constant with radius (outside of a "core" region, Rgc ^ ReS, that possibly has been dynamically depleted) if 
the recent surface-brightness profile of Caon et al. (1994) is used in the comparison rather than the older /is measurements 
of King (1978). The global specific frequency of NGC 1399 is similarly "average", Sn = 7.0 ± 0.6, if its GCS population 
is normalized to an absolute magnitude for the galaxy {Mv — —22.1) that is significantly brighter than the value more 
commonly adopted. The ratio Ed/Egtars is also constant with radius in this galaxy, and has a value equal to that found 
in M49. 
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M87, at the center of the Virgo Cluster, is much richer in gas than either of M49 or NGC 1399, and both Sn problems 
are in evidence there. A discrete, geometrical deprojection algorithm was developed to enable a direct comparison between 
the volume densities pd, Pstars, and pgas at all three-dimensional galactocentric radii in M87. It was found in §3.2 that 
Pel oc (pgas + Pstars)j bcyoud ^ 8 kpc, and that the constant of proportionality is the same as that in M49 and NGC 
1399. The cluster formation efficiency was constant with galactocentric radius inside each of these systems, and the same 
from one galaxy to another. Specifically, 

Cci ~ — = ^^^"l = 0.0026 ± 0.0005 . 

(Observed departures from this ratio in the central regions r^c ^ -Rcff of M87 and M49 may have resulted, at least in part, 
from the dynamical destruction of globular clusters there.) Many globulars in M87 are indeed associated with its X-ray 
gas, and this is in fact responsible for the appearance of both specific frequency problems. This single case leads to the 
basic expectation — which should be tested in as many individual systems as possible — that real discrepancies between Ed 
and Estars should appear only in gas-rich ellipticals that also have high global specific fequencies. 

The obvious generalization of these results is that most galaxies might have been subject to a single, common ed. 
This possibility was explored, and borne out, in §4. Observations of global Sn (i.e., total GCS populations and galaxy 
luminosities) in 97 giant ellipticals, brightest cluster galaxies, and dwarf ellipticals all were shown to be consistent with the 
predictions of a universal at the same level as that observed directly in M87, M49, and NGC 1399. The different scalings 
oi Sn with ii/,gai among normal giant ellipticals, central cluster galaxies, and dwarfs stem entirely from fundamentally 
different relations between iy,gai and in these different physical regimes. These relations have been obtained from 

scaling arguments that employ the optical fundamental plane of bright ellipticals; the observed correlation between X-ray 
and optical luminosities in early-type systems; and the fractional gas loss expected in a specific model (Dekel & Silk 1986) 
for winds from a single burst of star formation in dark-matter dominated dwarf ellipticals. Observations of A/tot vs. iy,gai 
still show some scatter about the mean prediction of constant ed, and it will be important to clarify the significance and 
the origins of this. However, there are no systematic deviations from a universal cluster formation efficiency, in galaxies 
that span four orders of magnitude in total luminosity and GCS population. 

Exactly the same ed was shown, in §4, to have obtained for globular clusters in the spheroid, or stellar halo, of the 
Milky Way; and it was argued that it appears also to apply now, to the formation of open clusters locally in the Galactic 
disk. This remarkable robustness, in which a single cluster formation efficiency was and is realized in a great variety of 
local and global environments, demands a very general theoretical explanation. From a more empirical standpoint, it 
bears directly on current questions involving GCS and galaxy formation and evolution on large scales. Some of these 
were discussed in §5. Of considerable interest is the fact that, with constant but Mgcs/Mstars not, the efficiency of 
unclustered star formation could not have been universal; field star formation must have been suppressed, in a sense, in 
dwarf ellipticals and at large galactocentric radii in the brightest, most gas- rich ellipticals. The possibility that these 
points could be understood in terms of feedback processes seems promising (§5), and should be explored in quantitative 
models. 

I am grateful to Bill Harris for providing the data plotted in Figures ^ and ^; to Ivan King and Chris McKee for helpful 
comments on an earlier draft of the paper; and to an anonymous referee for several valuable suggestions, including the 
addition of Figures and |l^. This work was supported by NASA through grant number HF-1097.01-97A awarded by the 
Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., 
for NASA under contract NAS5-26555. 
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DEPROJECTION ALGORITHM 

The surface density data in Tables ^, ^ and || lend themselves well to a geometrical deprojection of the sort that is 
often applied to the X-ray emission from ellipticals and galaxy clusters (e.g., Fabian et al. 1981; Kriss, Cioffi, & Canizares 
1983; Nulsen & Bohringer 1995). This technique makes only one assumption — that of circular symmetry — and is fully 
nonpar ametric. It also avoids the explicit differentiation of A'd that is required in the standard Abel-transform approach 
to deprojection. Such a differentiation can be performed if an ad hoc parametrization is adopted for the global form 
of the density profile from the outset, but it is clearly desireable to minimize any such a priori restrictions. Although 
sophisticated, nonparametric algorithms that resort finally to an Abel integral have also been developed (e.g., Merritt & 
Tremblay 1994), they present a greater degree of complexity than is necessary here; and, in any case, they seem poorly 
suited to density profiles that are constructed from multiple observations with different fields of view and statistical 
corrections for different levels of background contamination and photometric incompleteness. 

The assumption of spherical symmetry made here is not perfect (the isopleths of the M87 GCS, for one, are clearly 
non-circular in projection; see McLaughlin et al. 1994), but neither is it nonsensical: it takes azimuthally averaged surface 
densities and returns volume densities that are properly interpreted as angular averages at various galactocentric radii. 

A simple example of the basic method, which should serve to define its geometrical set-up and mathematical notation, 
is illustrated in Fig. |l^. The plane of the sky there is perpendicular to the plane of the page, and the two intersect in the 
fine z = Q. Thus, the fines on the right, labelled by projected radius i?o, and i?2, are fines of sight for an observer 
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Fig. 11. — Illustration of a geometrical deprojection algorithm. The plane of the sky is perpendicular to the page. Radii rg, ri, and r2 are 
three-dimensional quantities, and define spherical shells. Radii Rg, Ri, and R2 are projected quantities referring to cylindrical shells aligned 
along the line of sight; these correspond to circular annuli on the plane of the sky. 



situated off the bottom of the page. The bands Rq < R < Ri and Ri < R < R2, along with their mirror images on the 
left side of the Figure, are axial cross sections of cylinders that extend to z = ±00; in the plane of the sky, of course, 
the corresponding cross sections are circular annuli. The circular annuli shown here are not those on the plane of the 
sky, however; rather, they are equatorial cross sections of spherical shells with three-dimensional radii < r < ri and 
ri < r < r2- Evidently, these shells have been chosen such that vq = Rq, and so on. The origin r = (through which 
passes the axis R = 0) represents the center of a spherical distribution of globular clusters, which are supposed to have 
been counted up in discrete cylinders Rq < R < Ri and Ri < R < R2 on the sky, to form the surface densities Nci{Ro, Ri) 
and Nci{Ri, i?2)- From such data, it is possible to compute the fraction of clusters in each cylinder that actually reside in 
the spherical shells rp < r < ri and ri < r < r2, and thus to form the average volume densities nci(ro, ri) and nci(ri, r2). 

Suppose, for the moment (this will be relaxed just below), that r2 = R2 marks the outer edge of the GCS in question, 
i.e., nci{r > r2) = and Nci{R > R2) = 0. (It is assumed that the density Ni, of any uniform distribution of background 
objects has already been subtracted to form Nd from the raw data.) In this case, all of the objects seen in the outermost 



cylinder actually come from the outer spherical shell, i.e., from the hatched regions between i?i and R2 in Fig. 11, Since 
the number of globulars in this cylinder is just 

Ki{RuR2)^N,i{RuR2)xTT{Rl~Rl) , 

all that is required is a calculation of the volume of the shell ri < r < r2 which is intersected by the cylinder Ri < R < R2. 
Denoting this by t^nt(?'i, ^2] Ri, R2), the volume density is obviously 

, , Mci{RuR2) 
■nci{ri,r2) - 



Vint{ri,r2;Ri,R2) 



Moving inwards to the cylinder Ra < R < Ri, the total number of clusters observed there includes a contribution from 
the spherical shell ri < r < r2 (hatched regions between Rq and i?i), as well as one from the shell of interest, tq < r < ri 
(solid region). The volume density in the outer shell has already been determined, and the volume of intersection between 
that shell and the inner cylinder can be calculated [call this Vintifi, r2', Rq, Ri)], so the volume density in the inner shell 
is given by 

JVci{Ro,Ri) - nci{ri,r2)Vint{ri,r2; Rq, Ri) 

nci[rQ,ri) = — — — — , 

Vint(rQ,ri;Ro,Ri) 

where Afci{Ro,Ri) and Vint(rQ,ri; Rq, Ri) have the same meanings as the corresponding quantities in the first cylinder, 
Ri ^ R ^ i?2- 



Mclaughlin 



33 



Clearly, this procedure can be extended and applied iteratively, from the outside in, to the projected number counts 
in any series of concentric, non-overlapping, circular annuli; it is only required that there be no gaps between any two 
consecutive annuli. Thus, let the annuli be defined on the plane of the sky by projected radii Ro < Ri < R2 < ■ ■ ■ < Rm, 
i.e., set up m rings with Ro < R < Ri, R\ < R < R2] and so on. Given average surface densities A^ci in each of these 
rings, the volume densities rici in the corresponding spherical shells ro < ri < ?'2 < . . . < r,„ can be obtained once the 
volume of intersection between some generic shell Tj-i < r < rj and a cylinder Ri-i < R < Ri is known. A sketch similar 
to Fig. |ll| readily shows that this volume is, quite generally. 



]^nt(r'j_i,rj;i?j 



(Al) 



for J, J = 1,. 



This holds for any set of radii rj, Ri-i, Ri}, with the understanding that any term in parentheses 



which evaluates as negative is to be set to 0. An important application is to the case shown in Fig. [Ul, where, e.g., ri — Ri 
and r2 = R2 and the volume common to the shell Ri < r < R2 and the cylinder Ri < R < R2 isPint {^Rl t R2 T Rl: R2 ) — 
(47r/3)(i?2 — Rf)^^^- It is most appropriate to define shells such that r.; — Ri for all i, in which case the results given just 
above for m = 2 generalize to 
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Again, this assumes that any uniform background has already been subtracted to obtain the surface densities A'ci in every 
annulus i = I, . . . ,m. 

This equation is not quite final, however, since its derivation has been facilitated by treating the outermost annulus 
seen in projection, Rm-i < R < Rm (or the outermost spherical shell, Rm-i < r < R,n), as if it ended just at the 
"true edge" of some GCS. Of course, this will never actually be the case — if not because it makes little sense to associate 
sharp boundaries with such distributions in the first place, then because real observations and number counts can cover 
small fields of view within GCSs that are known to extend considerably further. This means that there is a second type 
of "background" contamination to be dealt with in any dataset, quite apart from that (already accounted for) due to a 
uniform distribution of unassociated galaxies and stars; corrections have to be made, in every annulus or shell, for the 
possible presence of bona fide globular clusters at radii r > Rm, beyond the last annulus for which a surface density can 
be directly measured. If the density distribution of such objects is denoted by nci(r), their contribution to the projected 
number density in a cylinder Ri^i < R < Ri is given by 



7?2 



RdR 



Ri- 



zq(R) 



n,i(r)dz- zo{R) = {Rl,-R^)^'^ 



(A3) 



which is just an average of the projected densities along every line of sight contained in the cylinder. [As usual, the 
coordinate z = (r^ — i?^)^/^ here measures distance along these lines of sight; z — Q defines the plane of the sky, as 
in Fig. 
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It is most convenient to express this "background" surface density as a fraction of the observed iVd in the 
outermost annulus Rm-i < R < Rm'- 
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-l,Ri) Rn 
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Ncl{Rm-l,Rm) R^ — R^ 



IrI.R'^R !Z(R) 
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(A4) 



where zq{R) is defined as in equation ( A_3), and the second equality relies on the basic definition of average surface density in 
a circular annulus of finite width. The advantage of this formulation is that the form of rtci(^) cannot be directly measured 
for r > Rm, so the calculation of TV^^'^^ necessarily depends on an extrapolation of some sort. Without any constraints on 

back,' E> . ^ Jl^^ ~> Nc\{Rm-l,Rm) 



this extrapolation, it is conceivable that an artificial and unphysical situation such as N]^y'^^{R. 



could arise. This and similar difficulties are avoided, however, by only ever evaluating N]^^"^^ in terms of / (which is 
guaranteed alwa ys to be <1 when i = m) and the known Nd at the largest projected radii observed. 

Equation (A4) should therefore be u sed to subtract N^^^'^^ from the total A'ci in every annulus before deprojecting the 
distribution according to equation ( A2). The corrected volume densities in the spherical shells i = 1, . . . ,m may then be 
written as 

3 

ndiRi-i, Ri) — n^^Ri-i, Ri) — — —j^ f{Ri-i, Ri)Nci{Rm-i, Rm) ■ (A5) 

Typically, when measured GCS surface densities extend to sufficiently large radius, an adequate estimate of / can be 
made by assuming a power-law distribution rid oc r~" in the outer reaches of the system, r > Rm', and in the case that a 
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is an integer, equation ( |A4| ) can be evaluated analytically. With a = 2, for example, 
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For a 
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And for a — A, 
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These three power-law exponents are expected to bracket the conditions that hold in real GCSs. 

Finally, the volume densities derived in this way are a series of averages of a continuous nci(r) in spherical shells of 
finite thickness. The question therefore arises as to how a single radius should be assigned to each shell so that the 
discrete distribution, nci(^i-i, Ri) vs. fi, represents the continuous one as faithfully as possible. (See Harris 1986 or King 
1988 for a discussion of this issue in the two-dimensional context.) Clearly, the answer is to choose each fi G [Ri-i,Ri] 
such that 



nc\{Ri-i,Ri) 



nci(r) dr = nc\{ri) 



(A9) 



R^-l 



If it is assumed that jid cx r^", the second equality is easily solved for r; in every shell. However, a will generally not 
be a constant from shell to shell; that is, a real CCS need not obey a single power-law density distribution over its entire 
extent (as is evident, for example, in Fig. || above). Moreover, a cannot be known a priori, and any attempts to derive 
it from the data must rely on fits to the average deprojected densities — fits which themselves depend on the values of 
r,;. Fortunately, a complicated iterative procedure can be avoided by noting that observed radial distributions suggest 
that <J a ^ 3 anywhere in a typical GCS. This in turn suggests, and numerical experiments confirm, that simply using 
a = 3/2 everywhere should lead to estimates that are acceptable substitutes for the exact solutions of equation (A£) 
in any realistic situation. Thus, the formula adopted here is 
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